Spectropolarimetric observations of Herbig Ae/Be Stars I: HiVIS 
spectropolarimetric calibration and reduction techniques. 
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ABSTRACT 

Using the HiVIS spectropolarimeter built for the Haleakala 3.7m AEOS telescope in Hawaii, 
we are collecting a large number of high precision spectropolarimetrc observations of stars. In 
order to precisely measure very small polarization changes, we have performed a number of 
polarization calibration techniques on the AEOS telescope and HiVIS spectrograph. We have 
extended our dedicated IDL reduction package and have performed some hardware upgrades to 
the instrument. We have also used the ESPaDOnS spectropolarimeter on CFHT to verify the 
HiVIS results with back-to-back observations of MWC 361 and HD163296. Comparision of this 
and other HiVIS data with stellar observations from the ISIS and WW spectropolarimeters in 
the literature further shows the usefulness of this instrument. 

Subject headings: Astronomical Techniques, Astronomical Instrumentation, Stars 



1. Introduction 

High-resolution linear spectropolarimetry mea- 
sures the change in linear polarization across 
a spectral line. Many models show this mea- 
surement is a useful probe of circumstellar en- 
vironments at small spatial scales. Circumstel- 
lar disks, rotationally distorted winds, magnetic 
fields, asymmetric radiation fields (optical pump- 
ing), and in general, any scattering asymmetry can 
produce a change in linear polarization across a 
spectral line such as Ha (cf. McLean 1979, Wood 
Brown & Fox 1993, Harries 2000, Ignace et al. 
2004, Vink et al. 2005a, Kuhn et al. 2007). These 
signatures can directly constrain the density and 
geometry of the circumstellar material and probe 
the near-star environment. Typical spectropolari- 
metric signals are small, often a few tenths of a 
percent change in polarization across a spectral 
line. Measuring these signals requires very high 
signal to noise observations and careful control of 
systematics to measure signals at the 0.1% level. 

This technique probes small spatial scales, typ- 
ically being sensitive to the region within 10 stel- 



lar radii. Even for the closest young stars (150pc), 
these spatial scales are smaller than 0.1 milliarc- 
seconds across and will not be imaged directly, 
even by 100m telescopes. Since the circumstellar 
material is involved in accretion, outflows, winds 
and disks, with many of these phenomena happen- 
ing simultaneously, spectropolarimetry can put 
unique constraints on the types of densities and 
geometries of the material involved in these pro- 
cesses. 

In order to address these issues, we built a spec- 
tropolarimeter for the HiVIS spectrograph on the 
3.7m AEOS telescope and performed a prelim- 
inary telescope polarization calibration (Thorn- 
ton 2002, Thornton et al. 2003, Harrington et 
al. 2006). Since then, we have made many im- 
provements. We have upgraded much of the hard- 
ware, calibrations, and reduction package. A large 
number of stellar observations have been partially 
reported in Harrington & Kuhn 2007 to illus- 
trate the capability of the spectropolarimeter. We 
present the reduction package, polarization cali- 
brations, and implications of these calibrations for 
stellar observations with this instrument. 



2. Improvements and A Dedicated Reduc- 
tion Package 

We have done several upgrades to the in- 
strument recently. The polarimeter has been 
remounted with better alignment and less vi- 
gnetting. A new tip-tilt guiding corrector has 
been installed and calibrated. A hard-mounted, 
adjustable dekkar and polarization calibration 
mount have been installed and allow for quick, 
easy, and repeatable calibration and adjustment. 
The failed CCD-ID20's were replaced with a sim- 
ilar copy, but an unresolved hardware/RFI issue 
caused this repair to fail. The original science 
camera failed and has been replaced by a new 
commercial Ik^ PIXIS camera and is working won- 
derfully. The calibration stage had lenses mounted 
that increased the flat-field flux by a factor of 6 
and the ThAr lamp flux by a factor of 150. These 
improvements greatly increase the ease and accu- 
racy of spectral calibration. 

The spectropolarimetric reduction pipeline can 
be broken into a few basic parts. 1) Combining 
flats, darks, biases. 2) setting up the coordinate 
grids as x-pixel, y-pixel, slit-tilt and wavelength 
calibration. 3) Calculating a flrst-pass spectrum 
by simply averaging all pixels at a given wave- 
length. 4) Calculating more spectra based on 
optimal-extraction routines, with cosmic ray rejec- 
tion and flux-weighting. 5) Doing the polarization 
calculations and 6) Simulating and quantifying the 
systematic errors. 

Combining flats and setting up the co- 
ordinate grid 

The first step is trivial - a clipped mean is per- 
formed on a number of flats to produce a flat 
field. Since we are doing high S/N measurements, 
we found that we needed a minimum of 20 flats 
at 10000 ADU to get good S/N. We decided to 
do the dark and bias calibration in a less-than- 
traditional way because of a slowly-fluctuating, 
non-repeatable background level, which is quite 
common with most detectors. We had an overscan 
region in the original science detector with the 
CCD-ID20's which allowed us to track the chang- 
ing dark levels throughout the night. The replace- 
ment Pixis 1024BR array was smaller than the fo- 
cal plane and is always exposed to light. Both the 
original Lincoln Labs CCD-ID20's and the Pixis 
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detector showed a significant (>few%) drift in the 
bias levels throughout a night. The Pixis cam- 
era had no overscan region, but the median val- 
ues of all the darks taken over many months did 
show variation in the overall background levels. 
These variations with this instrumental setup and 
background levels required a background subtrac- 
tion from each frame using the unexposed pixels 
outside the echelle orders as this directly tracked 
the background levels of the chip and reduced 
the systematic uncertainties. For our short ex- 
posures, the dark current was small, and is auto- 
matically corrected with the frame-by-frame back- 
ground subtraction. 

Since the instrumental setup and the order lo- 
cations on the CCD's has changed much over the 
last few years, and changed quite often during the 
instrument work done in the fall of 2006, setting 
up the coordinate grid requires more user-input 
than most dedicated reduction scripts. The user 
must supply a guess at the polarized order loca- 
tions, typically 5 points to fit a 3rd-order polyno- 
mial, and a bright stellar image taken that night 
to serve as a template for the curvature of the or- 
ders. Figure [T] shows a bright standard with the 
narrow Ha line near the center of the Ik^ array. 
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Fig. 1 . — The user-defined points serve as a start- 
ing point for the peak-finding routine that fits a 
3'"'' order polynomial to the shape of the order, 
as specified by a bright stellar template. The 
gaussian-fit centers are plotted as the jagged line. 
The polynomial fits to the peak intensity follow 
the stellar brightness peak quite nicely. 

The code sets up an initial polynomial fit to 
these user- input points, and then takes spatial 



(vertical) cuts at each spectral (horizontal) pixel. 
Each of these spatial cuts are then fit by a gaus- 
sian to determine a robust location of the peak 
intensity in the spatial direction across the CCD, 
reducing the influence of cosmic rays or chip de- 
fects. These peak intensity positions allow for a 
robust fit to the order shape on the CCD with a 
slope of around 5% (crossing a spatial pixel every 
18 or 20 spectral pixels). The orders are fit in- 
dependently and the shapes are the same to the 
0.1 pixel level (~0.2% of an order width). The 
order separation, caused by the Savart plate, is 
then calculated as the median difference between 
the two polynomial fits. For the Pixis camera, this 
was 58-59 pixels. The residuals between the cal- 
culated peak intensity and the polynomial fits are 
typically under half a pixel and show the expected 
oscillation with wavelength as spatial pixels are 
crossed. 

The star observed for the template is never 
completely centered in the slit from small guid- 
ing errors. Once we have these order shapes from 
the stellar observations, we then must adjust these 
polynomials to reflect the true slit center. This is 
done with the flat-field. First we create an aver- 
age spatial profile (slit illumination pattern) for 
the flat field, centered on the stellar template by 
simply averaging the flat field along the order tem- 
plate in the spectral direction. These spatial pro- 
files, shown in figure [2] clearly show the slit bound- 
aries, where the flat field intensity falls off sharply. 
The flat-field spatial cuts are cross-correlated to 
find the offsets from the stellar order-shape poly- 
nomial fits to the center of the flat-field spatial 
cuts. 

Once the offset from the stellar template to 
the flat-field region is calculated, the full coordi- 
nate grid can be computed with the order centers, 
edges, and background locations specified. Once 
the coordinate grid is set, we can then make a cal- 
ibrated flat fleld, adjusting for the blaze function, 
and wavelength dependence of the flat-field lamp 
(blackbody) . 

A flat-field is supposed to correct for pixel-to- 
pixel gain and fleld-dependent efficiencies only. 
The flat-field values should not introduce any 
blaze, lamp or spatial illumination dependencies. 
This means that the flat-fleld must be corrected 
for it's inherent spatial and spectral structure. To 
correct the flat fleld spectral dependence (horizon- 




Fig. 2. — The average spatial proflle of the two 
polarized flat-fleld orders across the Ha region. 
The width of the orders is set by by the dekkar 
width. The vertical lines mark the edge of the 
illuminated fleld we selected for use. The small 
tilt of the spatial cuts traces the inhomogenity of 
the illumination fleld at the slit. The difference 
in height between the two cuts corresponds to the 
intrinsic continuum polarization of the instrument 
downstream of the calibration stage. 

tal), a simple flat-fleld spectrum is made for each 
polarized order by averaging all pixels spatially at 
a given spectral position. This spectrum has a bit 
of intrinsic noise in it so a low-order polynomial flt 
is made. The flat-field image then has this spectral 
dependence removed. 

The spatial structure is corrected by dividing 
each spatial cut along the fiat-field spectrum by 
the average spatial profile, removing the spatial 
dependence of the fiat (inhomogeneous slit illu- 
mination) at each wavelength. This makes our 
fiat-field much closer to a true pixel-to-pixel gain 
measurement, shown in figure [3l 

2.2. Wavelength Calibration 

With the recent calibration-stage improvements 
we have implemented, it became possible to do 
very accurate wavelength calibration using low or- 
der polynomial fits to at least six clear, unblended 
lines, just over the Ha spectral region. With a 
>600x increase in ThAr flux (15x without the 
diffuser-bypass mirror), we can use the smallest 
slit (0.35") at highest resolution (R=49000) to 
clearly identify lines and compare the lines with 
standard line-lists obtained from the NOAO spec- 
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Fig. 3. — The calibrated flat field showing the cen- 
ter of each order and the width of each order. The 
region used for the background subtraction begins 
10 pixels above and below the order. The grey- 
scale runs from 0.95 to 1.05 showing the pixel-to- 
pixel structure. The region in the middle is ex- 
cluded from the reduction. 



The wavelength calibration is done in a two- 
step process. A Th-Ar spectrum taken with the 
0.35" slit is used to identify the lines. The Th- 
Ar spectra are calculated with the simple average 
over all illuminated pixels. The unblended lines 
are identified and the peaks are fit by a gaus- 
sian to determine the line centers. A 3'''^ order 
polyomial fit to the line centers is used to cre- 
ate a wavelength array. The NOAO atlas spec- 
trum, HiVIS spectrum, and a comparison AB- 
Auriga spectrum are shown in figure |4l The trian- 
gles mark the 6 lines used in the polynomial fits. 
Since all stellar observations were done with the 
1.5" slit, a 4-orientation sequence of ThAr spec- 
tra are taken with this slit and then reduced with 
the simple-average method. The resulting lower- 
resolution spectrum is cross-correlated with the 
higher-resolution spectra produced by the 0.35" 
slit to determine any shifts between slits and wave- 
plate orientations, the beam wander or wobble, 
which is a small fraction of a pixel. 



tral atlas. The new HiVIS detector has a wave- 
length scale of 0.037A per pixel (1.7 km/sec per 
pixel at Ha). Thus the rotation of the Earth pro- 
duces a i pixel shift and seasonal variations can 
approach 40 pixels. 
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Fig. 4.— The NOAO ThAr spectra is on the bot- 
tom. The calculated HiVIS 0.35" ThAr spectrum 
is in the middle. The triangle symbols mark the 
lines used in the polynomial fit. A sample unpro- 
cessed AB Auriga spectra is shown as the top line 
to illustrate the wavelength coverage of a single 
spectral order. 
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Fig. 5. — The spectral shift of the Th-Ar lines 
along the slit is calculated and used to set up a 
tilted coordinate grid. Main result is a 70% in- 
crease in resolution between reductions done with 
tilted and untilted coordinate grids (R=30000 to 
R=49000 for the 0.35" slit). 

The average slit tilt is then calculated to mod- 
ify the coordinate-grid for the reduction program. 
The tilt is a spectral-direction shift of ~7 pixels 
between the bottom and top of each order (across 
the 50 pixel slit width) . This shift is calculated as 
a cross-correlation of a spectral cut beween each 
spatial position along the order. This tilt is then 
implemented as a integer-pixel shift of the coor- 



dinate grid at each spatial point. There is a very 
significant increase in resolution when the tilted- 
coordinates are used. Figure [5] shows the spectral 
cuts before and after alignment, as well as the in- 
tensity of the line calculated with both tilted and 
un-tiltcd coordinate grids. The spectral resolu- 
tion derived from the Th-Ar line fits goes from 
R=30000 to R==49000 when using tilted coordi- 
nates. 

2.3. Calculating Rav^r Spectra 

Once the tilted coordinate grid, gain-table, and 
wavelength calibration is done, the data can be 
processed. Each exposure is divided by the gain- 
table (calibrated flat). A background spectrum 
is made by averaging the unexposed background 
region outside the polarized orders. This back- 
ground is noisy and so a linear fit to this back- 
ground is subtracted from each order to remove 
the dark current and bias of the detector. 

To create a simple, first-pass at the spectrum, 
every pixel at each wavelength was simply aver- 
aged to compute the intensity. At the same time, 
an average spatial profile (PSF) is constructed by 
averaging all these spatial cuts spectrally. While 
this method does not correct for cosmic rays, and 
does not weight the pixels by their fractional- 
stellar-flux or PSF, this simple average spectrum 
provides a starting point for more sophisticated 
reductions. 

2.4. Shift-n-Scale 

In many parts of the reduction package, we 
want to shift and scale one curve to another. We 
have developed a least-squares routine that allows 
us to calculate the shift and scaling coefficient be- 
tween two input curves. The routine assumes that 
the scale and shift from one curve to another can 
be written in terms of the input curve and it's 
derivative. To illustrate the routine we will show 
a shift-n-scale between the two polarized spectra, 
Ji shifted and scaled to /2, as below: 



all (x — b) 



(1) 



I is the intensity as a function of CCD pixel x, 
a is the relative gain between the two orthoganally 
polarized spectra, usually ^1, and b is the spectral 
shift in pixels, usually <0.8 pixels. The spectrum 
/i {x — b) can be expanded to first order as below: 



I2 = all — ab{dli / dx) 



(2) 



Rewriting the constants a = Ki and —ab = 
K2 and bringing all terms to one side gives the 
difference between the spectra I2 and /i scaled 
and shifted by a and b. The least squares routine 
minimizes the function below for all pixels (x): 



^{l2-Kih-K2dh/da 



(3) 
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Fig. 6. — An example of the shift-n-scale algo- 
rithm correcting for optical misalignments. A sin- 
gle spectrum, the bottom curve, was replicated 
and shifted spectrally by 1 pixel and then fed to 
the shift-n-scale routine. The top curve is the re- 
sulting large-amplitude (5% frj) fractional differ- 
ences. This curve follows the derivative and must 
be corrected. The shift-n-scale routine gave a scale 
of 0.999846 and a shift of 1.00742. The shift- 
n-scale interpolation causes some residual noise, 
since the input curves were identical and should 
have subtracted completely. This arises because 
of the smoothing done to calculate the derivative. 
The middle, mostly flat curve shows this residual 
noise doubled for clarity. The noise is at the level 
of 0.2% with some residual systematic error. 

The coefficients Ki and K2 are then used to 
calculate the scaled and shifted functions. An ex- 
ample of this routine applied to a spectrum of AB 
Aurigae is shown in figure [51 A single spectrum 
was replicated and shifted by one pixel and cor- 
rected by the shift-n-scale routine. The routine 
uses a smoothed profile to calculate the deriva- 
tive. This illustrates how a misalignment of the 
savart plate, simulated as a spectral shift of 1 pixel 



between the orthogonally polarized spectra, pro- 
duces a fractional difference of order 5% from this 
misalignment. This signature is a systematic error 
that follows the derivative of the spectral line. The 
scale and shift operation minimizes this derivative 
signal as shown (amplified) in the bottom line in 
the plot. The residual noise is at the level of 0.2% 
and has some residual systematic error in it from 
the smoothing used to calculate the derivative, but 
the derivative-based error has been attenuated by 
a factor of about 20. 

The shift-n-scale routine also has some nonlin- 
earity to it, and is best used for pixel shifts <0.8 
pixels. When shift-n-scaling the original curve 
to the simulated one, the routine reproduced the 
shifts to better than 0.05% for spectral shifts of 5 
pixels, far greater than the range needed. How- 
ever, when shift-n-scaling the simulated curves to 
the original curve, the calculated shifts had an er- 
ror of >0.1 pixels for spectral shifts >0.8 pixels. 
This sets the useful range of this routine. 
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Fig. 7. — This shows a cosmic-ray removal from a 
spatial cut. The cosmic-ray removal routine com- 
putes the deviation of the spatial cut from the av- 
erage line profile and identifies 4-a deviations in 
each spatial cut from the average spatial-profile 
(PSF) and replaces the cosmic ray hits with the 
value for the least-squares scaled average spatial 
profile. The smooth dashed curve is the computed 
least-squares scaled average spatial profile and the 
jagged curve is an individual spatial cut. Three 
pixels were found to be > 4(T away from the aver- 
age profile and were corrected. 



2.5. Optimal Spectral Extraction 

This first-pass spectrum and average spatial 
profile are projected back on the original spatial 
cuts with the shift-n-scale routine. The expected 
noise at each spatial pixel is found (using estimates 
of the gain and read- noise) . The standard devia- 
tion of each spatial pixel at each wavelength is 
then calculated and used to find any cosmic-ray 
hits as cosmic-rays (and chip defects) which will 
be strong devations (> 4tT) from the spatial pro- 
file. An example of this cosmic ray removal routine 
is shown in figure [7l The average spatial profile 
is the dashed line. The spatial cut at this wave- 
length shows a few pixels deviating strongly from 
the average spatial profile from a cosmic ray hit. 
These pixels are identified and replaced with the 
shifted-and-scaled spatial profile, greatly reducing 
the error in the spectrum. Once this cosmic ray 
removal routine is done, a spectrum is calculated 
as the simple average of all the pixels. The cor- 
rected spatial cuts are passed on to more optimal 
spectral extraction routines. 

The first of the optimal extraction routines uses 
the shift-and-scale routine to fit the average spa- 
tial profile (PSF) to each of the individual spa- 
tial cuts at each wavelength. This reduction pro- 
cess assumes that the spatial profile does not vary 
across the order, an assumption that all opti- 
mal reduction packages adopt (Home 1986, Marsh 
1989, Donati et al. 1997, Gushing et al. 2004). 
The calculated error in the shifts between the op- 
tional fit and the calculated order location is de- 
pendent on the flux, but is typically less than 0.1 
pixels. 

Another set of spectra are calculated as simple 
gaussian fits to the spatial cuts. While the PSF 
(spatial cut) is not entirely gaussian, and certain 
seeing conditions or tracking-errors lead to not- 
quite-gaussian spatial cuts (PSF's), the fits are 
usually robust, and match the shift-n-scale and 
simple-average specta well. The standard gauss- 
fit routine is applied to each spatial cut, with the 
initial estimates for the gaussian as: peak inten- 
sity is twice the computed average-spectra, cen- 
tered on the order, with a fwhm of one-quarter 
the order width and zero background. The rou- 
tine works well when the spatial cuts (PSF's) are 
roughly gaussian. 

Thus, we have three spectra with different 
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Fig. 8. — This shows a spectrum of AB Aurigae 
reduced with varying parameters to show the ro- 
bustness of the reduction program. Each curve 
represents a different reduction using order widths 
from 43 to 49, with or without flat fielding, and 
with or without tiltcd-coordinates. From top to 
bottom 1) tihed, 49 width, no flat 2) tilted, 43 
width, no flat 3) tilted, 43 width flattened 4) tilted 
45 width flattened 5) untilted 49 width no flat. 
The bottom curve shows the residuals from the 
unflattened untilted 49-width curve and the flat- 
tened tilted 43- width curve. 



systematic effects (simple-average, shift-n-scale, 
gaussian-fit) - which are all included in the po- 
larization analysis to probe our polarimetric re- 
duction codes and systematic errors inherent in 
these reductions. In practice, the exact choice of 
reduction parameters is unimportant. Figure [5] 
shows the simple-average method with and with- 
out a flat field or tilted coordinates, and varying 
order widths from 43 to 49 pixels (covering the 
edges of the illuminated region, sec figure [2]). The 
tilted coordinates do add a significant increase in 
resolution, but for the high signal-to-noise spec- 
tropolarimetric observations the spectra are re- 
binned to lower resolution anyways. The fiat field 
does increase the signal-to-noise of the intensity 
calculation, but as we show below, the basic po- 
larization calculations are independent of the flat 
field. 

The AEOS telescope is designed for military 



applications and does not have a standard siderial 
guiding telescope. There are significant tracking 
errors that can accumulate and there are offsets 
induced by manual telescope-operator corrections. 
The operation requirements set out by the AFRL 
require that we have no control over the guiding, 
and that their operators control the pointing of the 
telescope. There are a number of data set's where 
an operator's manual steps caused a non-gaussian 
spatial profile. Thus we calculate the spectra using 
all three methods, but the averaging and shift-n- 
scaling are most robust. 

These three different reductions vary slightly 
between each other. The gaussian-fit and simple- 
average methods vary by only a few parts-per- 
million when the spatial cuts (psf 's) are gaussian. 
However, the shift-n-scale method and the simple- 
average methods vary slightly across the line pro- 
file, giving an overall uncertainty in the shape 
of the intensity profile. However, the spectropo- 
larimetry is calculated as the difference between 
line profiles on the same image. Since the same 
reduction method is applied to each polarized or- 
der in the same way for each polarization measure- 
ment, it is only any differential error that would 
cause systematic errors in the spectropolarimetry. 
The reduction package uses all three calculation 
methods to assess the systematic errors that might 
be present in each method, or in the subsequent 
processing of the spectra. 

2.6. Continuum Subtraction 

Our focus is on spectral and spectropolarimet- 
ric variability. With the many instrumental polar- 
ization effects, we must do continuum-polarization 
subtraction to isolate the change across the line 
profile. Since spectral variability is always refer- 
enced to continuum, no spectrophotometric or flux 
calibration is done, though standards are observed 
for future work. 

The continuum is simply subtracted as a linear 
flt to the intensity on either side of the spectral 
line. This continuum subtraction normalizes the 
slope of the black-body of the star, removes the 
instrumental continuum polarization signals, and 
removes the instrumental blaze efficiency. Con- 
tinuum subtraction is done by averaging 20 pix- 
els sufficiently far enough away from chip edges. 
We found that there is a variable dark signal that 
curves the echelle orders near the edge of the chip. 



10-20 pixels wide, in a non-repeatable way, so this 
region is avoided. Cirrus, tracking errors, and 
wind-bounce all influence the throughput of the 
spectrograph and thus, even on a spectrophoto- 
metric night, a significant variation in the contin- 
uum flux is sometimes seen. 

The data taken before 10-27-2006 was imaged 
using an echelle order which placed the Ha line 
near the edge of the chip. This placement was 
done before we did some instrumental corrections, 
and before we remounted the Pixis camera to a 
different region of the focal plane. In the more 
recent data. Ha is closer to the blaze-efficiency 
peak and is much closer to the center of the focal 
plane. The continuum fits to the CCD-ID20 im- 
ages were not as accurate because the red wings of 
the Ha line were near the order edge, giving rise 
to a small (^1%) error in slope estimate. How- 
ever, this uncertainty is negligable for data taken 
on each night because the continuum is calculated 
in the same way for each instrumental setup on 
each night. This only produces an uncertainty in 
the night-to-night variability of the line profiles, 
and in the overall slope of the calculated contin- 
uum and continuum polarization. 

2.7. Possibility of Spectro-astrometry 

At this point it is useful to demonstrate the 
coordinate grid and basic properties of the reduc- 
tion package by demonstrating the wavelength de- 
pendence of the point-spread function properties. 
Several authors have argued that the change in 
centroid location (brightness peak) or full-width 
at half-max (fwhm) can be used as evidence for bi- 
narity in a set of spectroastrometric calculations 
(Beckers 1982, Bailey 1998, Baines et al. 2006). 
Spectroastrometry in principle is a technique that 
can detect unresolved binaries as a change in spa- 
tial profile width (psf) and location (centroid) 
across an emission or absorption line where the 
relative brightness of the two binary components 
significantly varies across the line. For example, 
if one star has an emission line and the other an 
absorption line, the centroid will shift towards the 
emitter in the core of the emission line. If the 
binary is significantly wide and bright, the fwhm 
would also decrease across the line as the relative 
contribution to the psf from the wide companion 
is reduced. This method makes use of the assump- 
tion that the psf across a line or spectral order does 



not change much, and that any significant change 
can be attributed to the source. HiVIS has been 
used to observe a number of young or emission-line 
stars as part of a large survey (Harrington & Kuhn 
2007) all of which have strong Ha lines. Given the 
reported detections, we will present the feasibility 
and capability of HiVIS to do similar work. 
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Fig. 9. — The spectroastrometry of MWC 166. 
The top box shows the center of a Gaussian fit 
to the PSF at every spectral pixel. The middle 
box shows the calculated FWHM for each spec- 
tral pixel. The bottom box shows the resulting 
normalized intensity. A systematic effect is seen 
across the line. 

Binaries are common among pre-main-sequence 
stars and this technique can easily be applied to 
every data set taken with this instrument. It also 
provides a good check on any systematic errors 
such as order mis-fitting or psf mis-fitting since 
any systematic error will show up with this mea- 
surement technique. Since the centroid of a spatial 
profile on most spectrographs can be calculated to 
a fraction of a pixel accurately in high signal-to- 
noise data, this technique has the capability to 
detect binarity to less than 30mas separations or 
a third of a pixel when our pixels are 130mas (cf. 
Bailey 1998, Baines et al. 2006). Normally, with 
a long-slit spectrograph, the slit of the instrument 
is rotated to a north-south line and an east-west 
line to fully sample the spatial changes. This al- 
lows one to establish the position-angle of the bi- 
nary system. For polarimetric reasons expanded 
below, we do not use the image rotator and keep 
this component fixed vertically during our obser- 



vations. This means that the position-angle ro- 
tates from exposure to exposure, but presenting 
the calculations is still useful to demonstrate the 
capability of the analysis software. 
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Fig. 10. — The spectroastrometry of 7 Cassiopeia. 
The top box shows the center of a Gaussian fit 
to the PSF at every spectral pixel. The middle 
box shows the calculated FWHM for each spec- 
tral pixel. The bottom box shows the resulting 
normalized intensity. No systematic effect is seen 
across the line. 

For each spectrum, the optimal spectral extrac- 
tion algorithms fit Gaussians and average spatial 
profiles (psf's) to the data at each wavelength, as 
was shown in figure [71 These spatial profiles can 
easily be compared with the nominal values for the 
average spatial profile's FWHM and centroid for 
each order. At each wavelength, a deviation for 
the FWHM and the centroid can be calculated. 
Because of small guiding errors, the centroid of 
each stellar observation does not coincide with the 
exact order center. Because of the tilt of the or- 
ders on the detector, this will cause the centroid 
of the spectrum to cross spatial pixels at differ- 
ent spectral pixel locations. In order to correct for 
this, our analysis package does a cross-correlation 
between the nominal order center location, as il- 
lustrated by the jig-saw shaped fit in figure [Tl and 
the centroid location of the actual observed spec- 
trum. Once this small offset has been corrected, 
the centroid location and the FWHM across the 
order is calculated. We will discuss some sample 
calculations from data taken the 20*'* of Septem- 
ber, 2007 below. 




Fig. 11. — The spatial profile or psf of P-Cygni 
on two separate exposures and three different re- 
ductions. One exposure was clipped (nearly sat- 
urated) to illustrate a systematic error in wave- 
length dependence of the psf. The top panel shows 
the psf at the emission peak. The triangles mark a 
gaussian fit to the non-clipped part of the psf. The 
bottom panel is the psf in the absorptive trough. 
The dashed lines in both panels show the Gaus- 
sian fits. Notice how the wings of the Gaussian 
are significantly above the noise in the low-flux 
absorptive trough. 

In Baines et al. 2006, the known binary MWC 
166 was reported to have a centroid shift of 
roughly 40mas EW and 15mas NW with corre- 
sponding FWHM changes of 40mas and 5mas re- 
spectively. This is tabulated in their paper as a 
total shift of 50mas at a PA of 287° where the 
known separation of the binary is 0.65" at a PA of 
298°. Figure [9] shows a HiVIS observation of this 
star. The observations had a peak of 1500 ADU 
and a continuum of 800-1300ADU across the order 
and a fwhm of 11.7 pixels (1.5"). The deviation 
we detect in FWHM and centroid offset are both 
around 40mas, both the correct magnitude for this 
system. 

We also present observations of 7 Cassiopeia 
in figure [10] at the same scale as MWC 166 to 
show a non-detection across a strong line. The 
peak intensity was near 26000ADU with contin- 
uum from 4000-6000ADU and a fwhm of 10.1 pix- 
els (1.3"). This non-detection shows that spurious 
signatures caused from the 1% Savart plate leak, 
non-linearity in the chip, spatial dependence of the 
focus or any other effects are not present across an 
emission line. 




Fig. 12. — The wavelength dependent psf prop- 
erties of P-Cygni for two separate exposures and 
three separate reductions. The three different re- 
ductions are tilted and untilted coordinates with 
width of 47 pixels, and untilted coordinates with a 
width of 41 pixels. The top box shows the center 
of a Gaussian fit to the PSF at every spectral pixel. 
The saturated exposure shows a smaller deviation 
than the unsaturated exposure for every type of re- 
duction meaning that more flux gives better fits. 
The tilted coordinates do not change the result. 
The middle box shows the calculated FWHM for 
each spectral pixel. The saturated exposure shows 
a SOmas increase in fwhm in the clipped part of 
the line. The unsaturated exposure shows a much 
greater 200mas decrease in fwhm in the absorptive 
component of the line. The tilted coordinates do 
not change the result, but using a smaller order 
width decreases the deviation, again suggesting 
that greater flux influences the fits. The bottom 
box shows the resulting normalized intensity. A 
systematic saturation effect is seen for all reduc- 
tions across the emission peak of line where the 
saturated exposure is underestimated. 



The star P-Cygni has a very strong Ha line. 
We observed this star with two different expo- 
sure times separated by 90 minutes. The reduc- 
tion was performed with three separate settings: 
a width of 47 for tilted and untilted coordinates, 
as well as a width of 41 un-tiltcd. This will be 
used to highlight changes in derived psf proper- 
ties when including tilts and neglecting low-flux 
psf wings. The altitude changed from 50° to 34° 
and the frame rotated by 28° between the ob- 
servations. The line has a normalized intensity 
of 0.194 in the absorption trough and 18.180 in 
the emission peak, or a flux ratio of 94. This 
huge range in flux and the difference in seeing 
serves to highlight systematics. One exposure is 
deliberately clipped (near saturated) to illustrate 
systematic errors that can occur with this tech- 
nique. The spatial profiles are shown in figure fTTl 
The peak intensity of the saturated exposure was 
44000ADU with continuum of 1800-3800ADU and 
only 400ADU in the absorptive component of the 
line. The calculated fwhm was 9.6 pixels (1.25"). 
However, a Gaussian fit to the unsaturated region 
shows that the true peak should have been closer 
to 48000ADU, or a roughly 10% clipping. The un- 
saturatd exposure taken an hour and a half later 
and at higher airmass had considerably less flux, 
16000ADU peak, continuum of 600-llOOADU, and 
only 130ADU in the absorptive component of the 
line. The calculated fwhm was 14.1 pixels (1.83"). 

The wavelength dependence of the psf proper- 
ties shows striking changes across the line, shown 
in figure [121 Both exposures and all reductions 
show a similar change in centroid location, sharply 
changing by about 50-150mas in the blue side of 
the absorptive component, then rapidly dropping 
across the center of the absorptive component to - 
20mas to -80mas in the red side of the absorption. 
The centroid then gradually rises over the emis- 
sion to a peak of 5-15mas just before the emission 
line returns to continuum on the red side of the 
line. The change is actually greater in magnitude 
for the unsaturated exposure with any reduction 
parameters. There are actually six curves in the 
top panel of figure [12] but the difference between 
order widths of 41 and 47 is completely negligible. 
The tilted coordinates show systematically lower 
deviations for both saturated and unsaturated ex- 
posures. In essence, more flux gives less deviation 
presumably from better fits. The tilted coordi- 
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nates more accurately reflects the true wavelength 
dependence of the orders and reduces the system- 
atic error. 

The fwhm also shows strong differences be- 
tween the different exposures and reduction pa- 
rameters. The unsaturated exposure with only 
140ADU peak flux shows a 200mas decrease in 
fwhm in the absorptive component of the line in 
both tilted and untilted coordinates with a width 
of 47 pixels. When the wings of the psf are ne- 
glected by using a width of 41, the calculated 
fwhm deviations are significantly smaller, near 
75mas, presumably because the wings of the psf 
have a noise contribution. This again shows that 
more flux leads to better fits and that large flux 
differences across a line lead to systematic effects. 
The saturated exposure shows an extra 50mas in- 
crease in fwhm across the emission peak in the 
clipped region of the line for all reductions. 




Fig. 13. — The wavelength dependence of the psf 
for T-Ori reduced with and without tilted coordi- 
nates. The top box shows the center of a Gaus- 
sian fit to the PSF at every spectral pixel. A 
systematic centroid shift effect is seen across the 
emission peak of line where there is extended, re- 
solved Ha emission. Using a tilted coordinate 
grid does not change the deviation. The middle 
box shows the calculated FWHM for each spec- 
tral pixel, smoothed for clarity. The tilted coor- 
dinates actually make the fwhm variation across 
the emission negligable, showing that this resolved 
emission does *not* contribute substantially to 
the fwhm. The bottom box shows the resulting 
normalized intensity. 

Another complication is extended emission. A 



change in fwhm and centroid can be caused by sig- 
nificant resolved emission around an object, such 
as our observations of T Ori. If there is no inspec- 
tion of the raw data, this emission may be falsely 
reported as a binary detection. Also - one cannot 
rule out the possibility of unresolved nebulosity as 
a significant contribution to the psf. The obser- 
vations of T Ori showed a peak signal of 600ADU 
with a continuum of 100-150ADU. The fwhm was 
calculated as 10.1 pixels (1.3"). The wavelength 
dependence of the psf for two reductions, using 
tilted and untilted coordinates, is shown in figure 
[131 The centroid shift is small, roughly 40mas and 
is independent of the order tilt. The fwhm changes 
by roughly ±80mas when using untilted coordi- 
nates, but is essentially zero when using tilted co- 
ordinates. This can arise when a vertical order ge- 
ometry intersects the fully-resolved emission near 
the edges of the order. When using properly tilted 
coordinates, the fully resolved emission becomes 
an added background that does not influence the 
Gaussian fitting. Figure [14] shows the raw expo- 
sure and the extended emission with the tilted co- 
ordinates shown as the slanted black line across 
the extended emission. 




Spectral PIxe 

Fig. 14. — The raw exposure of T Ori showing 
the spatially resolved extended Ha emission near 
the absorptive component of the stellar line. The 
lower beam has an untilted spectral slice shown, 
and the upper beam has a tilted coordinate grid 
overplotted. This illustrates a systematic error 
that may result from extended resolved nebular 
emission around a source causing a change in cen- 
troid and fwhm. 
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Further investigation of these effects is beyond 
the scope of this paper, but the analysis package is 
robust and subject to the same systematic errors 
as other analysis packages (Bailey 1998, Baines et 
al. 2006). In summary, this analysis package can 
reproduce results of others, specifically MWC 166. 
However, other systematic effects that always in- 
fluence spectroscopic data must be considered as 
well. This illustrates the robustness of our reduc- 
tion geometry and coordinate grid as well as the 
psf fitting and analysis by the optimal spectral ex- 
traction routines. 

2.8. Polarization calculation 

The simultaneous imaging of orthogonal polar- 
ization states also allows for a greater efficiency 
observing sequence and for reduction of atmo- 
spheric effects since polarization measurement can 
be done with a single image and both states are 
imaged simultaneously, with identical systemat- 
ics. Using the fractional polarizations measured 
as the difference between the orthogonally polar- 
ized spectra divided by the sum, q and u, we can 
calculate the percent polarization, P, and the po- 
sition angle, PA, as follows: 



O 1 ,a — b c — d ^ 1 , ,,, 

9=T = o(—^-:^-5)= 0(900 +9450) (4) 
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P = \J q^ + u^ 
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We wish to note that many people use the ratio 
method defined as follows: 
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The Ukirt infrared polarimetcr has yet another 
difference method published to calculate the po- 
larimetry on their IRPol2 website. 
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This method is a different version of the simple 
average method and can be reduced to a simple 
average where b and c switch places: 



Q l,a-c b-d 
1 2 a + c b + d 

When all three of these methods are applied to 
a data set, the difference between the calculated 
polarized spectra is entirely negligable. Since no 
method was significantly different than any other, 
we use the easiest to implement - a simple average 
applied to each individual exposure set. 

We wish to note that no instrumental polalar- 
ization ripple of the kind mentioned in Aitken & 
Hough 2001 or seen in ISIS (Harries et al. 1996) 
has been seen in our data. This ripple is thought 
to arise from Fabry-Perot type fringes in the wave- 
plate. Though our instrument uses a standard 
achromatic waveplate, a thin layer of polymer 
sandwiched between two glass substrates, no rip- 
ple near the 0.1% level has been seen. 

2.9. Aligning the spectra 

The Savart plate axis can never be perfectly 
aligned with the slit. This produces a small shift 
in the dispersion direction between the two polar- 
ized orders. This shift was calculated by cross- 
correlating the ThAr spectra from the two polar- 
ized orders. The derived shifts were 2 pixels in 
2004, 12 pixels in 2005 before the Savart plate re- 
installation, and 7 pixels after the 2006 remount- 
ing. 

To characterize any derivative noise, there are a 
number of ways we can align the spectra and cal- 
culate the polarization. The three ways are 1) no 
alignment 2) cross-correlation an a simple integer- 
pixel shift (which preserves noise properties) and 
3) integer-pixel shift followed by the shift-n-scale 
sub-pixel alignment. When no alignment is done, 
the spectra and polarization are simply calculated 
where the are on the detector. In the integer- 
pixel-shift alignment, we do a cross-correlation on 
the spectra to ensure that the spectra aligned to 
a single pixel, as a double-check on the wave- 
length calibration. Since any interpolation can 
effect and corrupt the noise statistics of the spec- 
tra, as shown by the shift-n-scale example in figure 
ini subpixel alignment must minimize the deriva- 
tive signal more than it introduces extra noise. 
To evaluate all three alignments, the polarization 
analysis is done on the three different spectra us- 
ing the three different alignment methods. These 
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nine polarization spectra are similar in their over- 
all structure. A sample of this is shown in figure 
1151 The polarization signature is around 1.5% and 
the differences between the different methods are 
very small. 



five error, as it is a derivative-of-derivative error 
(acceleration). We also calculated the sub-pixel 
shifts between the three different intensity calcu- 
lation methods using the shift-n-scale routine. We 
found the shifts to be a small fraction of a pixel. 
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Fig. 15. — Stokes q calculated 9 different ways and 
the differences between them- from the bottom, 
no alignment, integer-pixel alignment, subpixel- 
alignment, and the differences between all of them 
and the first one (simple average, no alignment). 
They're all quite similar in average shape. The 
changing S/N can be seen across the trough and 
peak of the PCygni line. The noise increases when 
any spectral alignment procedure is applied. 

The spectral shift (beam wobble) measured 
from Th-Ar lamp spectra taken at many waveplate 
orientations shows that there is a one pixel shift 
induced over a complete rotation of the waveplate 
with a roughly sinusoidal shape. The alignment of 
the waveplate with the savart plate causes a rota- 
tional offset between the rotation-stage zero point 
and the waveplate zero point of 85°. This means 
we operate our waveplate from rotation angles of 
85° to 152.5° near the crest of the sine-curve. 

Measurements of ThAr spectra taken each 
night, shown in figure 1161 show that the beam 
wobble is roughly 0.3 pixels between the four wave- 
plate orientations. The derivative signal caused by 
a 0.3 pixel shift for strong emission lines with con- 
trast of 5-20 is roughly 1-3%. It will be shown 
later that the derivative signal subtracts away in 
the dual-beam polarization analysis and that a 
beam wobble of 0.3 pixels produces a negligible 
error. This beam-wobble signal follows regions of 
curvature and is proportional to the initial deriva- 



Fig. 16. — The beam wobbles (shows a small wave- 
length shift) as the waveplate rotates. The wobble 
is calculated by applying the shift-n-scale routine 
independently to many lines in the Th-Ar lamp 
spectra. The wobble shows a sinusoidal depen- 
dence on waveplate angle with an amplitude of 
1 pixel. The 67.5° region of operation is near the 
crest of the sine-curve and hence has an amplitude 
of 0.4 pixels. 

The wavelength of each pixel for each horizon- 
tal order is assigned by a 3'''* order polynomial fit 
to thorium- argon calibration lamp images. This 
analysis is done independently on each order for 
the two orthogonal polarization states. Typically, 
reduction packages, such as Libre Esprit for the 
ESPaDOnS spectrograph on CFHT (Donati et al. 
1997) use the wavelength fits as the alignment 
method between the polarized orders. We have 
found that the simple average without alignment 
method works best because the flat-field cancels in 
the polarization calculation, leaving the measure- 
ment independent of the pixel-to-pixel gain. 

2.10. Rotating Frames 

Since HiVIS uses the alt-az coude focus, the 
spectropolarimcter's Q-U reference frame rotates 
on the sky as the telescope tracks. There are two 
main rotations: one with altitude, between the ter- 
tiary mirror and m4, as the altitude changes, the 
other with azimuth, between m6 and the coude 
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pickofF mirror. Normally, coude spectrographs use 
a 3-mirror image rotator somewhere downstream 
(in an pupil plane) to compensate for the tele- 
scope's motion and keep the orientation of the slit 
on the sky fixed. We have an image rotator, but 
since the rotator is 3 oblique reflections, we do not 
use this rotator to increase the accuracy of the po- 
larimetric calibration of the instrument. Changing 
the mirror orientations would change the polariza- 
tion induced by the mirrors. The orientation of 
the instrument on the sky is allowed to rotate and 
exposures are kept short so that the rotation is 
small. 

The analysis software takes the date and time 
of the exposure and calculates the altitude and az- 
imuth of the telescope during the exposure. Know- 
ing the orientation of the telescope mirrors, we 
can reconstruct the projection of the spectropo- 
larimeter's Q and U axes on the sky, and rotate 
all measurements to a common frame. A common 
convention is to set -l-Q to North on the sky. 

The rotation angle between the 4 individual ex- 
posures in each set is calculated as Position-angle 
-|- Altitude - Azimuth and is checked to quan- 
tify the inter-set smear. Very large rotation an- 
gles occured when an object is tracked near tran- 
sit (high azimuth rates). For measurements with 
small rotations, the polarization measurements are 
assumed to be accurate and are rotated to an ab- 
solute frame (of the first exposure) on the sky. 

For every data set, rotation in an individual 
exposure causes each exposure to become a mix 
of Stokes Q and U, reducing the efficiency of the 
measurement and adding some uncertainty to the 
Stokes parameters. Imagine an incident -l-Q beam. 
A rotation of the instrument of 22.5° during an 
exposure would cause the projected instrument 
frame to go from -|-Q to -|-U and the resulting mea- 
surement would only show 50% polarization, when 
100% was incident. 

Rotation between sets also induces systematic 
error. As the instrument rotates on the sky, dif- 
ferent amounts of time will be spent sampling the 
polarization states. As an extreme example, 22.5° 
rotation of the instrument between each individ- 
ual exposure would cause all images to be oriented 
on the sky with +Q whereas the reduction pack- 
age would assume the sequence was -f-Q, -|-U, -Q, 
-U. The systematic error with this high rotation 
rate causes complete loss of information about U 



while doubling the time spent measuring Q. It also 
causes the benefits of beam-swaping to be lost en- 
tirely and it causes the reduction methods to fail. 

2.11. Rebin By Flux 
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Fig. 17. — The rebin- by- flux routine adaptively 
averages pixels to equalize the S/N per spectral 
pixel. The number of physical pixels averaged per 
rebinned spectral pixel in the adaptive rebinning 
with a threshold of 4 is the bottom curve. In the 
P-Cygni absorption trough, 8 pixels are binned, 
but in the emission peak, there is no binning. The 
middle curve is the rebinned spectra which shows 
a highly compressed P-Cygni trough because of 
the high binning. The top curve is the flux per 
rebinned spectral pixel, with a minimum of 4 (the 
threshold) and a maximum of the emission peak. 

Since the expected spectropolarimetric signals 
are subtle differences in the polarized line profile 
shapes of less than one percent, and the few re- 
ported detections show polarization of at most 1%, 
the lines must be calculated with at least S/N> 
300 in all spectral-resolution-elements. 

Most young emission-line stars have large 
line/continuum ratio's (5-10). Some, such as AB 
Auriga, have PCygni profiles with emission-peak 
to absorption-trough ratio's of over 20. The star 
P-Cygni itself has a ratio of nearly 100. Even 
when saturating the peak of the emission line 
on the detector, a S/N of only 100 or so was 
achieved in the absorption trough of AB Auriga. 
A flux-dependent averaging procedure has been 
developed to bin pixels spectrally by flux, to en- 
sure that the S/N per spectral-pixel is roughly 
constant over the line at the expense of regular 
wavelength coverage. 
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Fig. 18. — The systematic error in a simulated 
data set caused only by varying noise. Random 
noise proportional to \/N was added to eight iden- 
tical copies of the Ha line of AB Aurigae, plotted 
above. The resulting Stokes parameters were cal- 
culated and show a noise of roughly 5% at con- 
tinuum and 2% in the emission line peak.. The 
raw Stokes q and u are the bottom curves, cen- 
tered about zero, with the noise obviously varying 
across the line profile. The spectropolarimetric 
signal, yg^ + i?, is the top curve and it varies 
by over 2% across the line profile. The variation 
is caused only by the varying noise, rising with 
higher noise. 



Another systematic error occurs when you cal- 
culate the degree of polarization for a spectrum, a 
squared quantity, with noise amplitudes that vary 
across the line. Since the average absolute val- 
ues of noisy data will be higher than clean data, 
and there will be significant noise variations across 
strong emission lines, varying noise can mascarade 
as a systematic spectropolarimetric error. Fig- 
ure [18] shows a simulated polarization measure- 
ment where flux-dependent noise was added to 
eight copies of a single line-profile. The noise 
varied from roughly 5% at continuum to roughly 
2% in the center of the emission line. The re- 
sulting Stokes q and u parameters were centered 
about zero, but with a varying noise amplitude. 
When these measurements are converted to a de- 
gree of polarization via ^J q^ -\- u^, this varying 
noise shows up as a systematic error in the degree 
of polarization. 



Fig. 19. — An example of the rebin-by-flux rou- 
tines applied to a data set. The noisy curves 
show Stokes q (top) and Stokes u (bottom) with 
a uniform spectral resolution, but an uneven S/N. 
The smoother curves near the middle have a more 
uniform noise per pixel and the wavelength de- 
pendence of the polarization is much more clearly 
seen. The average Ha line is overplotted. 

To rebin-by-flux, the average line profile for 
each data set (after alignment) is calculated. This 
provides the average flux and shape of all 8 lines. 
A user-defined threshold, typically 1-5, sets the 
minimum flux in each rebinned spectral pixel. The 
number of physical pixels to rebin at each wave- 
length is calculated and each individual spectrum 
is then rebinned accordingly. A rebinned wave- 
length array for each set is also calculated since the 
number of pixels binned in the spectral direction is 
different for each star and data-set. The polariza- 
tion analysis is then performed on these rebinned 
data sets as specified above and the pixel-to-pixel 
variation in the noise is greatly attenuated. Figure 
1 191 shows a comparision of raw and rebinned polar- 
ization spectra. The rebinned spectra only show a 
few resolution elements over the absorption trough 
with many covering the emission peak. The polar- 
ization signature in the absorptive part of the line 
stands out very clearly. 

3. Systematics - Alignment Error 

The calculated polarization is sensitive to sys- 
tematic errors in the alignment of the polarized 
spectra. There are two main alignments of the 
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spectra - between the polarized spectra in each in- 
dividual image, and between the specra from each 
individual image. 

The reduction code automatically makes derivative- 
error simulations for every data set to show the 
form of the systematic errors caused by misalign- 
ments. The misalignment between polarization 
states causes a derivative signal in each fractional 
polarization measurement u^d that can be at 
the few % level. As long as the shift between po- 
larization states is constant in each image, then 
the derivative signal is subtracted away when the 
second polarization measurement, :^~^: , is sub- 
tracted. The polarization state shift between im- 
ages is very stable. There is never more than a 
0.05 pixel shift difference between the calculated 
wavelength solutions of each polarized order as 
determined by ThAr line fits. 
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Fig. 20. — A beam- wobble polarization simula- 
tion - modeling the systematic error in polariza- 
tion caused by a spectral shift of the spectra in- 
duced by beam wobble or optical misalignments. 
As the waveplate rotates, the spectra of one orien- 
tation shift by a small fraction of a pixel to those 
of another orientation. This simulation greatly ex- 
aggerates the shift to show what a strong system- 
atic error would look like. The curves show the 
simulated Stokes Q spectra for beam wobble am- 
plitudes of 0, 0.5, 1.0, and 1.5 pixels. The largest 
shift shows a peak-to-peak amplitude of 0.4%. The 
systematic effect follows the curvature since beam 
wobble is a derivative of derivative. 

There is also a second order effect - beam wob- 
ble - that can cause a systematic polarization er- 
ror. If the spectra are misaligned between images. 



from the rotating waveplate or guiding errors in 
un-tilted coordinate systems, then a derivative- 
of-derivative effect can show up when the two 
polarization measurements are differenced. We 
simulate this effect in figure HOI by applying a 
orientation-dependent wavelength shift between 
the exposures that increases with waveplate an- 
gle. This simulates any beam-wobble introduced 
by misalignment of the waveplate. The induced er- 
ror traces the curvature of the emission line, with 
polarization errors showing up where curvature is 
greatest. Using ThAr spectra, we measured this 
wavelength shift to be '--^O.l pixels, depending on 
waveplate orientation, shown in 1161 This small 
wobble has a very small influence on the polariza- 
tion, as it takes a misalignment of several pixels 
to produce a noticable acceleration signal. 

The derivative simulations in this code show 
that the simple derivative error subtracts away 
completely if the polarized-order-separation is 
constant (which it is to 0.05 pixels). The code 
also shows that it takes a wobble of at least 2 full 
pixels between each orientation, for a total of a 
6-pixels misalignment to produce a 1% systematic 
error in a Stokes parameter for strong lines. Using 
our shift-n-scale routine, as well as the ThAr mea- 
surements, we have measured this wobble to be 
around 0.1 pixels and thus present only well be- 
low the 0.1% level proportional to the curvature 
(derivative of derivative) . 



Polarization Calibration 
Model 



The AEOS 



The next step in the instrument development 
was calibrating the instrumental polarization re- 
sponse. Absolute polarization at all wavelengths 
requires very careful calibration of the telescope. 
Internal optics can induce or reduce polarization 
(U,Q^I or I— >Q,U) in the beam, cause cross-talk 
between polarization states (QU^-V, V^QU), 
and also rotate the plane of polarization (Q^U, 
U^Q). For an instrument at the coude focus of 
an alt-az telescope, all of these effects are func- 
tions of wavelength, altitude and azimuth. For 
AEOS, there are five 45° reflections before the 
coude room, two of which change their relative 
orientation. This is particularly severe, and only 
a few put a spectropolarimeter after so many mir- 
rors. The calibration is difficult, but some have 
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tried this for many different telescope designs (cf. 
Sanchez Almeida et al. 1991, Kuhn et al, 1994, 
Giro et al. 2003) 

Since there are no spectropolarimetric stan- 
dards one usually observes polarimetric standards 
that have a precise value averaged over some band- 
pass such as Hsu & Breger 1982, Gil-Hutton & 
Benavidez 2003, Schmidt & Elston 1992, or Fos- 
sati et al. 2007. The calibration and creation of 
a telescope model is done by measuring unpolar- 
ized standard stars and polarized sources at many 
pointings to calibrate the effects of the moving 
mirrors. 

In general, the telescope will have different re- 
sponses to linearly polarized light since the absorp- 
tion coefficients of the optical components vary 
with wavelength and the differential absorption is 
a function of the angle of the incident polariza- 
tion with respect to the mirror/component axes. 
We used Zeemax models of the telescope to com- 
pute the assumed polarization properties of the 
coude focus given various optical constants for 
aluminum. These show the expected altitude- 
azimuth dependence of the mueller-matrix terms, 
though the amplitudes of the computed effects are 
significantly different than measured. While our 
models can calculate the full Mueller matrix as 
functions of altitude, azimuth, and wavelength, we 
cannot measure all of these terms. We can, how- 
ever, measure certain parts of the Mueller matrix 
given proper sources (polarized standards, unpo- 
larized standards, and scattered sunlight). 

Unpolarized standard stars are the first tool we 
used to calibrate the telescope. Observations of 
these standards at as many pointings as possible 
allows one to construct a model of the telescopes 
response to unpolarized light and to quantify the 
telescope-induced polarization (but not the depo- 
larization, cross-talk, or plane rotation which re- 
quires a linearly polarized standard). 

Polarized sources were the last calibration 
source we used. Using a polarized source allows 
one to measure the depolarization and rotation of 
the plane of polarization, all of which are, in prin- 
ciple, orientation dependent. Polarized standard 
stars are polarized in their continuum light, but 
the constant value quoted in standard references 
is an average polarization in a relatively line-free, 
wide bandpass image. The degree of polarization 
is usually not more than a few percent making in- 



strument calibration difficult with these sources. 
Another calibration method which we will explore 
below, is to use twilight (scattered sunlight) as a 
highly polarized source with a fairly well known 
degree and orientation of polarization. The polar- 
ization is certainly not constant and varies with 
atmospheric properties, time, and pointing, but it 
is a very bright and highly polarized source that 
has given us valuable insight into the polarization 
properties of the mirrors. 

4.1. Flat Field Polarization: Fixed Optics 
Polarization 

All of the moving mirrors in the AEOS tele- 
scope are located upstream of the polarizing op- 
tics (analyzer) . As an illustration of the polariza- 
tion induced by the reffections between the halo- 
gen flat field lamp {Tefj ~ 2900i4r) and the spec- 
tropolarimeter, figure [2T] shows the polarization 
reduction applied to a set of flat field frames with 
the 4 waveplate orientations taken in 2004. The 
flat field lamp is located just after the coude room 
entrance port, immediately after the coude pickoff 
mirror (m7). The lamp is reflected off a diffuser 
screen and is assumed to be unpolarized. 

The flat fields show a significant polarization 
which changes with wavelength. At 6500A, the in- 
duced polarization is dominated by the -t-U term, 
but the U and -Q terms become equal at 90. 
The position angle varies almost linearly across all 
wavelengths while the degree of polarization rises 
from 2% at 650nm to 6% at 950nm. Since the 
remounting of the optics and a movement of the 
image-rotator, the values of the fiat field polariza- 
tion have changed, but this illustrates the magni- 
tude and wavelength dependence of polarization 
induced by the 11 fixed reflections in the optics 
room. 

We also have tested the polarization of the 
common-fore-optics. A linear polarizer was 
mounted (with a mask) at various points along 
the optical path, from the calibration stage just 
after the coude port to just in front of the slit. 
A linear polarizer was measured to have 98.2% 
polarization at the slit, 98.3% polarization after 
the k-cell, 97.4% before the k-cell, and 97.7% after 
the flat field screen. This shows that the linear 
polarization does not vary much between the be- 
ginning of the spectrograph optics and the slit. 
Another polarizer mounted at the slit was mea- 
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Fig. 21. — Flat field polarization analysis, aver- 
aged 200:1 for clarity. Stokes q is shown as dia- 
monds , u is shown as triangles, and the degree of 
polarization is the top curve, shown as -I- symbols. 
The squares represent the position-angle of polar- 
ization (PA) with the y-axis scale on the right. 
The polarization (top curve) begins at 650nm be- 
ing dominated by u (triangles), but begins to rise 
as -q (diamons) gets larger. The PA rotates by 
about 25°, going from nearly 55° or mostly +u to 
being nearly 70° or a mix of -q and +u near (-3,4). 

sured to have 98.4%. This polarizer was put at 
the calibration-stage and was measured to have 
98.1% 97.7% 97.8% and 98.0% polarization for 
angles -03.2, 44.0, 87.5, and 134.2 degrees in the 
slit's QU frame when manually rotated to 0, 45, 
90, and 135 degrees. The measured angles match 
completely within the manual uncertainty. This 
shows that the degree of linear polarization is pre- 
served through the spectrograph optics. 

The flat field also shows a difference in beam 
intensities (vignetting or beam transmission dif- 
ference), in addition to it's intrinsic polarization. 
The polarization in a single frame, computed as 
^ shows +q at -4.55% -q at 1.36% +u at -3.63% 
and -u at 0.45%. These values are not symmet- 
ric about zero, as would be the case if there was 
only intrinsic polarization. To make the polar- 
ization values symmetric about zero, a constant 
difference of 1.59% between the two beams is ap- 
plied. This was calculated as 1.588 and 1.591 for q 
and u respectively as the null polarization spectra 
(sometimes called the check spectra: -|-q -I- -q). 
Restated, this means the ±q ±u spectra are sym- 
metric about 1.59%, ie that the rotation of the 



waveplate swaps polarization around a constant 
difference in beam intensity of 1.59%. Subtract- 
ing the 1.59% difference in beam intensities from 
all single-frame polarizations ^^ makes these po- 
larization spectra symmetric about 0% refiecting 
the true polarization of the beams. 

We note that the old fiat field polarization mea- 
surements showed a polarization of 2.78% at a PA 
of 55.0° (-0.93% 2.57% qu) for 65671. This is not 
what we measured with the new setup. This is 
expected since the image rotator was moved, the 
fiat-field lenses were installed, and the new dekkar 
is slightly smaller than the old one, changing the 
illumination pattern. The new fiat field polariza- 
tion is just as repeatable and robust as under the 
old configuration. 

Table 1: Flat field polarizations 



Name: 


+q 


-q 


+u 


-u 


Pol %: 


-4.546 


1.364 


-3.629 


0.452 


Pol-1.59%: 


-2.956 


2.954 


-2.039 


2.042 



The calculated flat fleld polarization after back- 
ground correction is -2.95% for q and -2.0% for u 
calculated with the average method as 0.5* ( -f-q - 
-q ) giving a degree of polarization of 3.59% at an 
angle of 287.3 degrees. We performed this test on 
a sequence of 40 flat field frames and each individ- 
ual exposure set gives the same answer to 0.01% 

We did a test of the dichroic for the LWIS 
mounted inbetween m4 and m5. The dichroic is 
roughly 12" square and causes an estimated 12% 
reduction in throughput. We measured the twi- 
light polarization without, with, and again with- 
out the dichroic mounted. The insertion and 
removal takes less than one minute so the twi- 
light polarization and PA can be assumed con- 
stant throughout the ^5 minute observation. The 
measured polarization values were nearly identical 
with and without the dichroic. 

4.2. Unpolarized Standard Stars - Tele- 
scope Mirror Induced Polarization 

Observations of unpolarized standard stars al- 
low us to measure the variation of the telescope- 
induced polarization as a function of pointing. 
Since the incident light is unpolarized, any ob- 
served polarization signature directly traces the 
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Fig. 22. — The measured polarization of unpolar- 
ized standard stars for the Ha region pfotted in the 
QU plane. The average value is -2.74% q -2.38% u 
or an average of 3.8% at 295°, which is essentially 
the flat-field polarization. There was no obvious 
correlation with these points and either altitude 
or azimuth. 

polarization induced by the telescope - the IQ lU 
mueller matrix terms. Many unpolarized stan- 
dards have been observed in November 2004, June 
& July 2005, Sept 2006 - Jan 2007, and June 
2007 totaling 224 data sets. Since the camera 
was changed in September 2006 and a better wave- 
length solution was obtained, only the 126 obser- 
vations after September 2006 were used. The cal- 
culated degree of polarization for each unpolarized 
standard observed is averaged to a single measure- 
ment over the entire Ha region. Note that all our 
unpolarized standards showed no significant line 
polarization, and had a statistical error of less that 
0.05%, typically 0.02%. These polarization mea- 
surements are shown on the q-u plane in figure [22l 
The measurements are centered around the flat- 
field value of -2.95%q -2.04%u. When adjusted 
for the flat field polarization, they show a range of 
polarizations from 0.5% to 3.0% 

The observations can be plotted on an altitude- 
azimuth plane to make an all-sky induced-polarization 
map. The altitude-azimuth coverage, shown in fig- 
ure [231 was fairly uniform from altitudes of 20 to 
80 degrees. There is a slight gap low in the North 
that comes from a lack of unpolarized stars in 
this region - the highest declination star we are 
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Fig. 23. — The altitude-azimuth coverage for all 
the unpolarized standard stars used in making the 
telescope-induced polarization maps. 




Fig. 24. — The measured polarization of unpolar- 
ized standard stars for the Ha region in the alt-az 
plane. The range of polarization is 0.5% to 3% for 
130 independent observations. 



aware of is at +50 (HR7469). These observations 
are interpolated to a regular alt-az grid. Figure 
[Ml shows the all-sky map of the telescope- induced 
polarization. This surface does not show a signifi- 
cant alt-az structure. To understand the apparent 
variability in the map, we did a test of an un- 
polarized standard star, HR8430, as a function 
of the secondary mirror focus. There is a signifi- 
cant variability, or order 1% as the focus changes 
shown in figure [25l This magnitude of variability 
with secondary focus, combined with the induced- 
polarization map suggest that the telescope in- 
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Secondary Focus 

Fig. 25. — The measured polarization of an unpo- 
larized standard star - HR8430 - for the Ha region 
as a function of the secondary mirror focus. There 
is a strong dependence on focus value. The best 
focus value was 5 and the width of the psf more 
than doubled for focus values of 1 and 10. 

duced polarization map has a very small alt-az 
dependence, of the same order as the focus depen- 
dence. 

4.3. Scattered Sunlight - A Linearly Po- 
larized Source for Calibration 

Scattered sunlight is a strong linearly polarized 
source and is reasonably described by a simple 
single-scattering Rayleigh model (Coulson 1988). 
Singly-scattered light is polarized orthogonal to 
the scattering plane with the degree of polarization 
proportional to sin^O. The twilight polarization is 
highest at a scattering angle of 90° which reaches a 
maximum 90° from the sun on the North-Zenith- 
South great circle when the sun is in the West. 
The degree of polarization in a Rayleigh atmo- 
sphere can be simply described with this equation: 



d — d„ 



sing 



1 + cos g'^ 



cos g = sin tz sin t cos p + cos t^ cos t 



(10) 



(11) 



Where g is the angular distance between the 
pointing and the sun, t^ is the solar elevation, t is 
the pointing-to-zcnith angle, and p is the pointing- 
to-solar meridian distance. Contained the variable 
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Fig. 26. — The Rayleigh model - the degree of po- 
larization of scattered sunlight with the sun set- 
ting in the west projected onto the alt-azi plane. 
The polarization is highest where the image is 
white. The maximum polarization was set to 80% 



g is the altitude-azimuth dependence for a given 
location of the sun. At sunset, t^ ^0, simplifying 
the model. Our observations were done with the 
sun lower than about 15°, typically from 45 min- 
utes before sunset to just after sunset. Using this 
simplification, we can rewrite cosg as sin i sin p 
and simplify the degree of polarization as: 



d = d„ 



sm g 



1 + (sin f sin p)^ 



(12) 



This gives us a simple square-sinuosoid in an- 
gular distance from the sun (g) modulated by the 
pointing-to-zcnith and meridian distances. The 
North-Zenith-South meridian is 90° from the sun 
and represents the maximum polarization d„iax- 
The plane of polarization would always be paral- 
lel to the altitude axis (vertical). In the absence 
of optical depth effects (which depolarize the light 
near the horizon), the polarization would always 
be djnax on this meridian. On the meridian from 
the East-Zenith- West, the degree of polarization 
would vary as sin^, being zero in the East and 
West and dmax at the zenith. The plane of polar- 
ization would always be perpindicular to the EZW 
meridian and aligned with the azimuth axis of the 
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telescope (horizontal) . Figure [26] shows this sim- 
ple Rayleigh model for the twilight polarization 
projected onto the altitude-azimuth plane. 

At our high-elevation (10000ft) the dry, clear 
atmospheric conditions are particularly good for 
using a Rayleigh model because of the low optical- 
depth and high degree of polarization. Near sunset 
at normal telescope pointings, scattered sunlight 
has a known position angle and a degree of polar- 
ization of typically 30-80%. Observing scattered 
sunlight at many different altitudes and azimuths 
allows us to check the instrument's response to lin- 
early polarized light as a function of mirror angles. 
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Fig. 27. — Scattered sunlight polarization at an 
elevation of 75° for azimuths of NSEW from top 
to bottom. Polarization of scattered sunlight is 
maximum at 90° scattering angle and at sunset, 
this is the arc from North to South through the 
Zenith. The polarization spectra have been av- 
eraged 200:1 for ease of plotting. The symbols 
are: q=0, u=A, P=^/Q^ + U^=+. P is the top 
curve with Stokes q and u are the more sinuosoidal 
curves with a strong rotation of the plane of po- 
larization (itan~^-). 

This scattered light is not a perfect source how- 
ever. The degree of polarization is a strong func- 
tion of the scattering properties of the air. In par- 
ticular, aerosol content (salt, dust, water, etc), op- 
tical depth (multiple scattering changing with al- 
titude and horizon distance) and reflections from 
the Earth's surface (land and ocean) can signif- 
icantly change the observed linear polarization 
(Lee 1998, Liu & Voss 1997, Suhai & Horvath 



2004, Cronin et al. 2005). Our site is surrounded 
by a reflecting ocean surface, and it is unknown 
what effect this has on the twilight polarization 
at our site. Since the sun was up during most 
of our observations, the polarized light reflected 
off the ocean surface has potential to complicate 
our measurements. However, we expect this ef- 
fect to be small because the maximum polarization 
at zenith obsered over many nights was relatively 
constant, being ^ 75%, while the cloud covering 
the ocean below changed significantly. The broad 
trends in the degree of polarization we measured 
are in agreement with the singly-scattered sky po- 
larization model, but the linear polarization ob- 
served by many researchers has been shown to be 
functions of time, atmospheric opacity, aerosol size 
distributions, and composition, all of which vary 
significantly from day to day (Coulson 1988, Liu & 
Voss 1997, Suhai & Horvath 2004). We have no ac- 
curate, quantitative simultaneous observations of 
the twilight polarization to compare with our ob- 
servations, but there are still many useful things 
to learn from measurements of a highly polarized 
source. 
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Fig. 28. — The calculated angle of polariza- 
tion {^tan~-^ ^) for scattered sunlight at altitudes 
[30,50,65,75] and azimuths [N,E,S,W] with the sun 
. The 52 observations were taken June 26th to 
July 7th 2005 with somewhat irregular sampling 
on this alt-az grid. 

We previously reported observations of twilight 
at many different pointings from June 26th to July 
7th 2005 in Harrington ct al. 2006. We obtained 
52 twilight polarization measurements (416 spec- 
tra) on an altitudes- azimuth grid of [30°, 50°, 65°, 
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75°, 90°] and [0°, 90°, 180°, 270°] covering 6500A 
to 9500A with the old science camera. Due to time 
constraints and a much longer read-time, we were 
unable to obtain more than 5 measurements on 
any one evening, and there is significant variability 
at single pointings of over 5% night-to-night. We 
then compared these measurements to expected 
sky polarizations (Coulson 1988, Horvath et al. 
2002, Pomozi et al. 2001). 

These spectra showed varied wavelength depen- 
dences of polarization and position-angle. A sam- 
ple of polarization spectra in the north, south, 
east, and west at 75° elevation are shown in fig- 
ure [27l They have been binned to 200 times lower 
spectral resolution, five points per spectral order, 
for ease of plotting. The figure shows the north- 
south and cast-west pairs are almost identical even 
though the north and west spectra were taken one 
night after the south and east spectra. Most of the 
52 spectra showed stronger polarization changes in 
degree and angle at shorter wavelengths with the 
degree of polarization becoming more constant at 
longer wavelengths. However, the rotation of the 
plane of polarization was a very strong function of 
wavelength and pointing, with sometimes 90° of 
rotation from 6500 A 9500 A). The strong rotation 
of the plane of polarization for all 2005 twilight 
measurements is shown in figure 1281 These curves 
do not sample the alt-az grid evenly, but they all 
show strong, relatively smooth rotation with wave- 
length. Some show a very strong dependence be- 
tween 6500A and 7000A. Even if the degree of 
polarization becomes more constant, the rotation 
of the plane of polarization continues with wave- 
length. It is hard to explain this as anything but 
the telescope. 

We then did much more extensive twilight ob- 
servations with the new hardware in June 20th- 
22nd 2007. The AEOS dome walls were raised, but 
the aperture was open, allowing us to observe only 
altitudes above 60°. We obtained 109 polariza- 
tion measurements over the evenings of June 20th, 
21st, and 22nd, typically spending an hour to com- 
pletely map the sky. On all three nights we cov- 
ered at least one complete map and we always had 
multiple overlapping observations. Since the sun's 
azimuth at sunset was 296° during these three 
nights, we aligned our azimuth coverage to re- 
flect this as our west zero-point. The observations 
spanned 4:46-5:19UT on 6-20, 4:11 to 5:01UT on 



6-21 and 4:23 to 5:01UT on 6-22. Over this time 
frame, the sun's altitude dropped by about 10° 
and the azimuth changed by 5° . The night of the 
20th, we obtained one complete grid at altitudes of 
60,75, and 89 degrees elevation at azimuths of 26, 
116, 206, and 296 degrees, the equivalent of NESW 
in the simple rayleigh model. On the subsequent 
two nights, we extended our azimuth coverage to 
include azimuths 71, 161, 251, and 341, or the sec- 
ondary coordinates (NE, SW, etc). We obtained 
nearly three full patterns on the 21st and nearly 
two full patterns on the 22nd. 




Fig. 29. — The average measured polarization of 
scattered sunlight as a function of altitude and 
azimuth during sunset for all observations on the 
20th-22nd. The twilight polarization is assumed 
to follow a simple rayleigh-scattering model with 
an azimuth-independent maximum at zenith. The 
saw-tooth form at high altitudes is entirely the 
telescope. 

As with the unpolarized standards, the scat- 
tered sunlight polarization measurements can be 
plotted on the altitude-azimuth plane. Figure [^ 
shows the measured degree of polarization pro- 
jected on the sky for the June 2007 measurements 
around the Ha line. There is a striking 4-ridge 
structure in the degree of polarization that im- 
mediately stands out at the zenith. The oscilla- 
tions are strongest at high altitude and decrease 
in amplitude as the altitude falls, becoming double 
peaked by an altitude of 60°. The position- angle 
also shows strong variation with pointing. 

Figure [3D] shows all the measurements plotted 
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Fig. 30. — The measured polarization of scattered 
sunlight as a function of altitude during sunset. 
The polarization should rise with altitude, but 
one rising trend (along cardinal NESW azimuths) 
and one falling trend (along intermediate azimuths 
SW, SE, NW, NE) is distinctly seen. The dashed 
line represents the region all the points would oc- 
cupy for a simple Rayleigh model with a maximum 
degree of polarization of 75%. 



versus altitude. The polarization at the zenith is 
at its maximum in the simple rayleigh model and 
should not vary at all with changing azimuth. The 
spread in the values at each azimuth shows the ef- 
fect of the telescope on polarized hght. On all 
three nights there were observations at a pointing 
repeated an hour after each other. The change 
in polarization in that hour was as much as 10% 
(rising by 10% at zenith as the sun set). However, 
the change in polarization from the telescope was 
roughly 5 times this amount in azimuth. In an 
ideal setting, the points would all converge to a 
maximum polarization (~75%) at the zenith while 
maintaining a spread from low to high at lower al- 
titudes. The dashed line in the figure shows the 
simple Rayleigh model range for a maximum po- 
larization of 75%. This is obviously not the case, 
and two things stand out. The first is that there 
are no points at maximum polarization at lower al- 
titudes. In the simple rayleigh model, observations 
to the north and south should be just as strongly 
polarized as the zenith. The second is that the po- 
larizations measured at the zenith range from 10% 
to 80%. Since the measurements at the zenith rep- 



resent a 2° patch of sky, and the time variability 
is of order 10%, the telescope must cause an ef- 
fective depolarization from 75% to 15% for some 
pointings. 
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Fig. 31. — The measured polarization of scattered 
sunlight at the zenith during sunset as the tele- 
scope rotates in azimuth. The first measurement 
is at 292° at 4:16UT, wrapping around and fin- 
ishing at 316° at 4:56UT. The polarization should 
be constant with all azimuths, but a sinosoidal 
variation is seen. This signature is caused by the 
crossing and uncrossing of the m6-coude pickoff 
mirror pairs. The minima at 160° and 340° are 
around 10% whereas the minima at 70° and 250° 
are around 30%. The maxima are always 65-75% 
with a systematic, nearly linear increase with time. 

In order to further examine these surfaces, we 
performed more twilight measurements at the 
zenith on July 27th and 28th 2007. Both nights 
were relatively clear but with slightly higher hu- 
midity. On the 27th, we pointed the telescope 
at the zenith and spun the azimuth axis over a 
number of points, covering 360° in 12° steps with 
three extra points overlapping. Figure [31] shows 
the computed degree of polarization and figure [32] 
shows the position angle of polarization in the in- 
strument frame. The first measurement is at 292° 
at 4:16UT, wrapping around and finishing at 316° 
at 4:56UT. Since the telescope never deviated from 
the zenith, the sinuosoidal oscillation in the degree 
of polarization is entirely caused by the telescope. 
In the hour it took to make these measurements, 
the sun had dropped several degrees in altitude. 
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increasing the overall degree of polarization at the 
zenith. This is causing the points in the right side 
of the figure to not close the sine-curve. The plot 
of position angle has been folded back on itself 
by 180° since the Stokes parameters are inher- 
ently 180° ambiguous. The position angle should 
go through two complete 90° rotations. The po- 
sition angle does this rotation, however it does 
not do it linearly! The deviation is at least 30° 
from linear and the deviation is nearly identical 
over the hour of the measurement, in contrast to 
the 10% change seen in the degree of polarization 
over the same time period. Since the rotation of 
the plane of polarization matches so well over the 
entire hour of observations, we are confident that 
this rotation is entirely the telescope. Also, dur- 
ing the hour this experiment was performed, the 
overall intensity of the sunlight decreased linearly 
by 50% with no dependence on azimuth, show- 
ing that there is minimal absorption of different 
polarization states. 
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Fig. 32. — The measured position- angle of polar- 
ization of scattered sunlight at the zenith during 
sunset in the instrument frame as the telescope 
rotates in azimuth. Since there is a 180° ambigu- 
ity in the Stokes parameters, the data have been 
folded by 180° with the measured angles match- 
ing almost perfectly. The polarization should be 
rotate linearly with azimuth, but an extra varia- 
tion of up to 30°is seen. This signature is caused 
by the crossing and uncrossing of the m6-coude 
pickoff mirror pairs. The relative azimuth is set to 
160° to match the minimum polarization seen in 
the twilight measurements. 



Another thing to notice is that the minimum 
polarization at 160° and 340° is close to 10%, 
whereas the minimum at 70° and 250° is more 
like 30%. On the 28th, we did a more in-depth 
study of the polarization minima in order to more 
fully quantify the minima. Wc observed azimuths 
of 138-180° and 232-274° in 3° steps spending 13 
minutes at each group at 4:26UT and 4:40UT re- 
spectively. Figure [33] shows the degree of polar- 
ization and measured instrumental position angle 
for these observations. The azimuth's have been 
shifted by 90° to center them. The deep minimum 
shows a very strong PA flip through the minimum. 
The shallow minimum shows a much weaker angle 
dependence. 
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Fig. 33. — The measured degree of polarization 
and position-angle of polarization of scattered sun- 
light at the zenith during sunset near the max- 
imum depolarization azimuths. The solid lines 
show the degree of polarization with the scale on 
the left. The dashed lines show the instrumental 
position angle of polarization with the scale on the 
right. The polarization measurements at azimuths 
of 232-274°, shown in triangles, have been shifted 
by 90° to fit on this plot. 

In summary, the twilight observations show 
that the telescope can cause a 75% linearly po- 
larized source to depolarize to 15%. The telescope 
can also rotate the plane of polarization by at least 
30° with pointing. The telescope also causes the 
rotation of the plane of polarization with wave- 
length by up to 90° from 6500A to 9500l. All 
of these effects are strong functions of pointing. 
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The 2005 measurements suggest that the depo- 
larizing effect greatly diminishes at longer wave- 
lengths, but that the rotation of the plane of polar- 
ization continues. These are very strong telescope 
polarization effects that must be considered when 
discussing results with this spectropolarimeter. 

4.4. Zemax Model of Telescope Polariza- 
tion - Mueller Matricies 

We constructed a model of the telescope's po- 
larization response using Zemax to compare the 
predicted polarization effects with our observa- 
tions. Zemax software allows one to propagate 
any number of arbitrarily polarized rays of light 
through any optical design. Each material or sur- 
face can be given any desired material property 
(thickness and complex index of refraction). We 
wrote programs to compute the Mueller matrix of 
the telescope for any altitude or azimuth by track- 
ing polarized light through our telescope with a 
range of mirror angles. We traced the polariza- 
tion properties of the telescope from the pupil to 
the first coude focus, just after the coude pickoff 
mirror (mirror 7). This allows us to get a quali- 
tative idea of the altitude-azimuth dependence of 
the polarization effects induced by the relative ro- 
tation of the altitude and azimuth mirror pairs. 

We started by using the optical constants for 
aluminum from the Handbook of Optical Con- 
stants of Sohds II (Pahk 1991) and ended up test- 
ing a range of refractive indicies from 1.5-2.5 for 
the real part and 6-10 for the imaginary part. 
These values are near the range of values mea- 
sured for aluminum in the green to near infrared 
regions. For instance, in the Handbook of optics 
the index goes from n-t-ik of (0.77, 6.1) to (1.47, 
7.8) to (2.80, 8.5) for wavelengths of 500nm to 
650nnm to 800nm respectively (Bass 1995). Also, 
Giro et al. (2003) find an index of (0.64, 5.0) in 
a V filter (530nm) using a model of the Nasmyth 
focus of the alt-azimuth TNG telescope to fit ob- 
servations of unpolarized standard stars. This is 
entirely within the reasonable range. 

To compute the Mueller matrix, we input six 
pure polarization states ±Q, ±U, and ±V at 400 
pupil points each, propagate the light through the 
telescope at a given pointing, and average the re- 
sulting polarized light at the image plane. Each 
reflection in the Zeemax models are given the same 
complex index of refraction. The six results enable 




Fig. 34. — The Zeemax calculated induced po- 
larization, calculated as 100*(IQ^+IU^)2 at the 
coude focus for an index of refraction n=1.5-8i. 
This shows the expected continuum polarization 
is of order 2% to 15% as a function of pointing. 
A higher complex index reduces the degree of po- 
larization, but does not significantly change the 
pointing dependence - all indicies give a double- 
sine with azimuth, rising with altitude. 



you to calculate the 16-element Mueller matrix for 
the telescope at a given pointing as the transfer 
from each input polarization to each output po- 
larization. These six calculations were computed 
for each pointing in 1° steps in altitude and az- 
imuth mirror angles to project the Mueller matrix 
elements on to an all-sky map. These maps were 
computed for a number of refractive indicies to de- 
termine the dependence of the telescope polariza- 
tion on the assumed real and complex index com- 
ponents. From these sky maps the expected am- 
plitude and pointing-dependence of the telescopes 
polarization effects can be seen. For notational 
simplicity, we will assume a Mueller matrix nota- 
tion as in this equation: 



M = 



/ // QI UI VI \ 

IQ QQ UQ VQ 

lU QU UU VU 

\ IV QV UV VV J 



(13) 



The complex component of the refractive in- 
dex was the dominant term - the higher the com- 
plex index, the shallower the wave can penetrate 
the aluminum and hence the less significant the 
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Fig. 35. — The Zeemax QQ Mueller matrix term 
at the coude focus for an index of refraction n=1.5- 
8i. The matrix element has been multiplied by 100 
for axis clarity. This term shows how the overall 
reference frame of the instrument rotates with az- 
imuth - a double-sine curve with azimuth is ex- 
pected for a telescope with no polarization effects. 
The slight bow seen at mid-altitudes reflects rota- 
tion of the plane of polarization by the telescope, 
causing small deviations from a perfect double- 
sine curve. 



has two angles where the mirrors are aligned (par- 
allel and antiparallel), we expect a double-ridge 
structure in azimuth. 
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Fig. 36. — The Zeemax calculated rotation angle 
in degrees at the coude focus for an index of re- 
fraction n=1.5-8i. The angle is calculated as the 
difference between the ideal telescope (with com- 
plex index -inf) and the simulated telescope. The 
angle is calculated as half the inverse tangent of 
QQ on QU. 



polarization effects. We also ran a model with a 
complex index of 100, meaning essentially no alu- 
minum penetration, and the models predicted per- 
fect polarization response by the telescope. The 
real component of the index had a much smaller 
effect, with higher indicies meaning more perfect 
polarization response, as expected. 

The telescope induced polarization from un- 
polarized incident light was computed as the IQ 
and lU Mueller matrix terms added in quadrature, 
100*(IQ2-|-IU2)5 , shown in figure[34l The induced 
polarization is expected to be minimal when the 
mirrors are crossed and maximal when they are 
aligned. The rotation of the plane of polarization 
can be calculated as the deviation of the frame 
rotation from the ideal frame rotation. The ideal 
frame rotation is shown in figure [35l This rotation 
angle is computed as ^atan{QQ/QU) with the 
large complex refractive index. Figure[3nishows an 
example of the telescopes rotation of the plane of 
polarization. The rotation angles should be maxi- 
mal at intermediate angles. Since the azimuth axis 



These models proved ultimately to be quan- 
titatively inaccurate. Table [5] shows the range 
of Mueller matrix elements for a range of opti- 
cal constants as well as the range of polarization- 
plane rotation angles and induced polarizations. 
The predicted telescope-induced polarization (IQ 
and lU) were much greater, 2-28% depending on 
pointing and index, than the measured 0.5%-3.0% 
polarization in unpolarized standards. Even the 
double-sine shaped pointing dependence predicted 
by these induced-polarization models is masked by 
other effects. 

The rotation of the plane of polarization was 
predicted to be quite small, 3° for a complex in- 
dex of 10 to 8° for a complex index of 6 with an 
altitude-azimuth dependence illustrated in figure 
1361 However, we observed very strong rotations, 
certainly at least 30°, as previously seen in figures 
[28l and [32l The model predicted the depolariza- 
tion terms, QI UI and VI terms to be nearly iden- 
tical to the induced polarization terms, IQ, lU and 
IV. This predicted moderate to low depolarization, 
which clearly contradicts the twilight observations 
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in figure [31] 

The rotation of the plane of polarization was 
also predicted to be only mildly dependent on 
wavelength since the complex index rises with 
wavelength from 6 at 500nm to near 11 at lOOOnm 
to over 20 at 2000nm. The deviation of the QQ, 
UU, QU and UQ terms from an ideal rotation ma- 
trix was small, at most 8° with a complex index of 
6, and the rotation decreases with a rising complex 
index. The indicies of refraction for aluminum do 
not change by much over this wavelength range. 
In the Handbook of optics, the index of refrac- 
tion goes from (1.4, 7.8) at 6500l to (2.8, 8.5) 
at 8OOOI and then to (1.2, 11.2) at IISOOA (Bass 
1995). The complex index of refraction is always 8 
or higher, meaning the model will predict less than 
4° rotation of the plane of polarization. This also 
has been measured to be completely false by the 
twilight observations since the plane of polariza- 
tion rotates by over 90° for some pointings, shown 
previously in figure [28l 



Table 2: Zeemax Mueller Matrix Terms 



rvlueller 


Index 


Index 


Index 


Index 


Term 


(1.5, 6) 


(1.5, 8) 


(1.5, 10) 


(2.5, 8) 


IQ, Qi 


-0.17 to 0.28 


-0.10 to 0.17 


-0.06 to 0.11 


-0.16 to 0.26 


lU. UI 


-0.22 to 0.22 


-0.13 to 0.13 


-0.09 to 0.09 


-0.20 to 0.20 


IV, VI 


-0.04 to 0.08 


-0.02 to 0.03 


-0.01 to 0.02 


-0.03 to 0.05 


VQ, QV 


-0.22 to 0.77 


-0.17 to 0.63 


-0.14 to 0.53 


-0.16 to 0.60 


VU, uv 


-0.87 to 0.62 


-0.75 to 0.50 


-0.64 to 0.41 


-0.71 to 0.47 


vv 


0.41 to 1.00 


0.64 to 1.00 


0.76 to 1.00 


0.66 to 1.00 


Rot (°) 


-8.1 to 7.6 


-4.4 to 4.3 


-2.8 to 2.7 


-4.2 to 3.7 


% Pol 


5.8 to 28.2 


3.4-16.4 


2.2 to 10.1 


5.2 to 26.5 



We are unable to measure circular polarization 
with our current instrument configuration, but the 
circular-induced term (IV and VI) were predicted 
to be very small, less than 0.1 for the lowest com- 
plex index and less than 0.02 for the highest. The 
circular cross-talk terms (VQ, VU, UV, QV) were 
all predicted to be very significant, running up to 
nearly 1 as a function of pointing and index. This 
also implies that significant V sensitivity will be 
lost at certain pointings, as illustrated by the VV 
term being less than 50% for some pointings in 
figure EH 

There is also further evidence that the index 
of refraction is significantly different for this tele- 
scope at longer wavelengths. Circular cross talk 
has been measured from other instruments using 
the AEOS telescope. The Lyot project uses a near- 
infrared (JHK) coronographic imaging polarime- 



ter mounted in a different coude room (Oppen- 
heimer et al. 2003). This room changes the zero- 
point of the azimuth mirror, but does not change 
the optical path. They also utilize the adaptive- 
optics system for the telescope, which adds sev- 
eral reflections just before the coude pickoff and 
outputs a coUimated beam instead of a converg- 
ing beam. The instrument first has to change 
the beam diameter by a number of oblique folds. 
The polarization analysis is then performed in the 
coUimated beam just after the coronograph with 
liquid-crystal waveplates and a WoUaston prism. 
Despite the differences between our optics and 
wavelengths, they also found very significant Q, U, 
V crosstalk that was pointing dependent in H band 
(1600nm) images (Oppenheimer et al. 2008). The 
total polarization \J (f ^ u^ ^ v^ for AB Aurigae 
was constant over a range of pointings though the 
individual terms varied. Since the complex index 
is predicted to be over 16 at this wavelength the 
Zeemax models would produce very small cross- 
talk even with several added mirrors. 

The aluminuum oxide coating that forms on ev- 
ery mirror is a major source of uncertainty. Since 
we are unable to model the effect of aluminum ox- 
ide on the mirror surfaces, which always form, are 
significant, and possibly variable in time, the mod- 
els are not expected to be quantitative. The hand- 
book of optics states that the oxide layers typically 
reduce the optical constants by 25% in the IR, 10- 
15% in the visible, and not much in the UV (Bass 
1995). This reduction is small compared to the 
range of input parameters used for this analysis, 
however, with so many mirrors, we suspect the 
contribution to be significant. Furthermore, our 
model does not include any of the spectrographs 
common fore optics. 

The polarization measurements of these com- 
mon optics described above showed that the de- 
gree of polarization remains nearly identical to 
the input polarization and that there are no 
orientation-dependent effects. Complete linear 
polarization at the coude focus (via a linear po- 
larizer at the calibration stage) is completely re- 
produced by the polarimeter. This suggests that 
these common foreoptics are not that significant in 
their overall effect on the degree of polarization for 
any input angle. However the measurement does 
not constrain an orientation-independent rotation 
(a PA offset). The flat field measurements shows 
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3.6% polarization, showing that the polarization 
induced by the common fore optics is also small, 
even though there are 8 reflections including the 
three highly-oblique image-rotators reflections be- 
tween the calibration stage and the slit. This also 
suggests that most of the severe polarization ef- 
fects are caused by the telescope mirrors. These 
Zeemax models are only an illustrative guide of 
what polarization the relative rotation of the mir- 
rors can induce and what the pointing dependence 
may look like. 




Fig. 37.— The Zeemax calculated VV Mueller 
matrix term at the coude focus for an index of 
refraction n=1.5-8i. The matrix element has been 
multiplied by 100 for axis clarity. This term shows 
a signiflcant loss of circular polarization (V) to 
other polarization states at certain pointings. 

We had claimed in Harrington et al. 2006 that 
the Zeemax model gives justiflcation for neglect- 
ing circular polarization. The circular polariza- 
tion terms IV and VI were predicted to be quite 
small, typically less than 0.05. The VQ, QV, UV 
and VU terms, which describe incident circular 
polarization becoming measured linear polariza- 
tion were much larger, typically 0-0.7. Since most 
astronomical sources have nearly zero circular po- 
larization and the V signatures across Ha in our 
stars has been measured to be very small we con- 
cluded that neglecting circular polarization was 
justifled (cf. Catala et al. 2007 or Wade et al. 
2007). We have now shown these models to be 
quantitatively inaccurate and are uncertain of the 
telescope-induced circular polarization. However, 



even though these effects are expected to be se- 
vere, since we are observing line polarization, even 
severe circular crosstalk will not impart more than 
a continuum-offset to our polarization measure- 
ments since they have negligable wavelength de- 
pendence across a single spectral line. The linear 
polarization detections we report are typically an 
order of magnitude greater than the circular signa- 
tures reported elsewhere (cf. Harrington & Kuhn 
2007 and Wade et al. 2007) and these cannot be 
the result of cross-talk. 

4.5. Implications for spectropolarimetry 

The telescope polarization effects can be quite 
severe and complicated. The effect on stellar spec- 
tropolarimetry is signiflcant, but useful observa- 
tions can still be performed with this instrument. 
The telescope can not induce any polarization ef- 
fects across a spectral line because the wavelength 
dependence of the induced polarization is quite 
weak at these wavelength scales. This means 
that any observed spectropolarimetric signature at 
some wavelength in a spectral line must be real, 
and at least at the observed magnitude. The mea- 
sured signature could have been depolarized, ro- 
tated, translated, or be a cross-talk signature from 
the stars circular polarization, but the wavelength 
of the change in the line cannot be affected and 
the amplitude serves as a lower limit on the true 
spectropolarimetric effect. 

The effect of telescope polarization is best il- 
lustrated in the so-called qu-loops. Polarization 
spectra are typically presented as either q and u 
spectra in iqu plots, or as degree of polarization 
and position angle spectra in ipt plots. Some au- 
thors augment these with plots of the polariza- 
tion in qu-space, ie plotting q vs u (cf. Vink 
et al. 2005b or Mottram et al, 2007). When 
a polarization spectrum is plotted in qu space, 
where each wavelength becomes a point in q and 
u, the telescope's polarization effects can be sim- 
ply represented as rotations, translations, shears, 
and amplitude- attenuations. Figure [38] shows an 
example of a q and u spectrum of AB Aurigae 
plotted in the qu plane. The spectra have been 
smoothed heavily to remove noise and clarify the 
example. The overall intensity is plotted in the 
top left panel. The q and u spectra, shown in 
the top right and bottom left respectively, show 
strong spectropolarimetric effects in the absorp- 
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Fig. 38. — A sample intensity, q, and u spectra 
of AB Aurigae, averaged heavily and plotted with 
wavelength and in the qu plane. The top left box 
shows the intensity. The top right and lower left 
show q and u respectively. The bottom right shows 
a plot of q vs u, illustrating a qu-loop. The contin- 
uum polarization has been set to zero, so a cluster 
of points representing this continuum for a knot 
around zero. The loop is caused by the change of 
q and u across the absorptive component of the 
line. 

tive part of the line profile. In the qu plane, when 
plotting q vs u, these spectra show a loop. The 
polarization of the continuum outside the line cre- 
ates the so-called continuum knot at zero. Since 
the continuum polarization has been removed from 
these spectra, each point on the graph from out- 
side the spectral line plots near zero forming a 
cluster of points. The spectra show a change of 
-1.4% and 1.7% for q and u respectively, which 
occurs near the red side of the absorption trough. 
These changes send the qu-loop off to at a PA of 
135° to a maximum amplitude point (-1.4, 1.7). 

Figure[35]shows an example of a telescope depo- 
larization effect in qu-space. A qu loop is shown 
with an initial amplitude of about 0.4%. With 
increasing telescope depolarizations, the telescope 
reduces the polarization of the source and the am- 
plitude of a qu loop decreases. Telescope-induced 
continuum polarization would cause these loops to 
cluster about some non-zero point. Rotation of the 
plane of polarization could rotate this loop to an 
arbitrary angle. In effect, the telescope can shift, 
reduce, or rotate a qu-loop, but it cannot create 
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Fig. 39. — A simulation showing how the telescope 
nullifies a qu-loop as it becomes an increasingly 
strong polarizer. The telescope attenuates polar- 
ization at some direction until one polarization 
state is completely extinguished. If the qu-loop 
is aligned with the position-angle of the polarizing 
telescope, the qu-loop is completel nullified. How- 
ever, to extinguish a qu-loop, the telescope must 
be nearly 100% polarizing. 



one. 

4.6. Calibration Summary 

We have illustrated the polarization effects of a 
altitude azimuth telescope at the coude focus. The 
moving mirrors can induce, reduce, and rotate po- 
larization as functions of position, pointing, and 
wavelength. We measured these effects in a num- 
ber of ways. We used linear polarizers mounted at 
various points in the optical path to polarize the 
light and measure the output. We used unpolar- 
ized standard stars as a method to measure the 
telescope-induced polarization. We used twilight, 
a highly polarized source, to check the depolariza- 
tion and rotation of the plane of polarization. 

Our measurements with a linear polarizer 
mounted behind the spectropolarimeter showed 
that we can detect linear polarization downstream 
of the spectropolarimeter with an efficiency better 
than 95% in the useful wavelength range of the 
linear polarizer. The polarizer mounted just after 
the coude pickoff mirror, in the calibration stage at 
the entrance of the room showed that the common 
fore optics produce a negligable depolarization for 
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all linear polarization orientations. 

To do a complete calibration of the induced po- 
larization, measurements of unpolarized sources 
have been made at many pointings. The unpolar- 
ized star measurements have been projected back 
onto the sky so that a map of the telescope re- 
sponse has been made for the induced polariza- 
tion to allow us to calibrate other sources of inter- 
est. The unpolarized standard star measurements 
showed that the telescope induces polarization of 
order 0.5%-3.0% for the Ha region in the new con- 
figuration. 

The scattered sunlight degree of polarization 
measurements showed very strong effects of the 
telescope on a highly polarized source. The mea- 
surements at the zenith showed that the telescope 
can cause an effective depolarization, from 75% 
to 15% for the Ha region. This effective depolar- 
ization was seen to be much less severe at longer 
wavelengths in our 2005 observations, but we have 
not repeated those measurements with the new 
configuration. The plane of polarization also ro- 
tates quite strongly. The wavelength dependence 
of the rotation of the plane of polarization was also 
shown to be strong across all wavelengths, even as 
the depolarization effect decreased significantly. 

The telescopes polarization properties are quite 
severe and comphcated. The effect on spectropo- 
larimetry is significant, but the telescope can not 
induce any effects across a spectral line. In effect, 
the telescope can shift, reduce, or rotate a qu-loop, 
but it cannot create one. Thus, any spectropolari- 
metric signature we detect is a useful lower limit 
on the magnitude of the polarization change across 
a line, and a good indicator of the wavelength 
range of the polarization effects. 

5. Instrument comparison: ESPaDOnS 
vs. HiVIS 

As a verification and extension of our HiViS 
observations, we observed two targets with HiVIS 
and with another spectropolarimeter, ESPaDOnS, 
mounted on Mauna Kea's 3.6m Canada-France- 
Hawaii telescope (CFHT). ESPaDOnS is a fiber- 
fed spectropolarimeter with a resolution of 68000, 
roughly 5 times higher than HiVIS in it's typi- 
cal lowest-resolution polarimetric mode (Donati et 
al. 1999, Manset & Donati 2003). The spectropo- 
larimeter has a completely different design with 



the polarizing optics mounted after the Cassegrain 
focus in a triplet-lens collimated beam. The 
instrument uses a dual-beam design with three 
fresnel-rhomb retarders and a woUaston prism 
which produces the two orthogonally polarized 
beams. The two polarized beams are then fo- 
cused onto two fibers (with two additional for sky 
that are unused in polarimetric mode). The typ- 
ical exposure sequence is 8 exposures at differ- 
ent Frenel-rhomb orientations, 4 each for Q and 
U. The telescope is equatorial and thus does not 
have any frame-rotation issues as the fibers are 
always at a fixed orientation on the sky. This 
makes the instrument very useful for long integra- 
tions. The fiber entrances do however introduce 
a strong instrumental continuum polarization, re- 
quiring a polynomial subtraction similar to HiVIS, 
and making continuum polarization studies impos- 
sible with this instrument. This instrument has a 
very different setup from HiVIS and as such is a 
useful cross-check. 
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Fig. 40. — The spectropolarimerty of MWC361 
from HiVIS on the nights of June 18-21st 2007 
is shown on the left and from ESPaDOnS on 
the nights of June 23 & 24th 2007 on the right. 
The MWC361 spectropolarimctry from HiVIS has 
been arbitrarily rotated by-eye to align the HiVIS 
qu-loops with those from ESPaDOnS. Both data 
sets have their continuum polarization removed in 
the respective dedicated reduction scripts by low- 
order polynomial fits. Both have been rebinned to 
lower spectral resolution for clarity. 

We observed MWC361 with the ESPaDOnS 
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Fig. 41.— The spectropolarimerty of HD163296 
from HiVIS on the nights of June 18-21st 2007 
is shown on the left and from ESPaDOnS on the 
nights of June 23 &24th 2007 on the right. No 
rotation has been apphed to the HiVIS data. The 
HiVIS data spans a range of pointings and ex- 
posure times. There is a resemblance between 
some HiVIS and ESPaDOnS observations, but sig- 
nificant telescope polarization effects cause the 
amplitude of the HiVIS signature to be smaller. 
Both data sets have their continuum polariza- 
tion removed in the respective dedicated reduc- 
tion scripts by low-order polynomial fits. Both 
have been rebinned to lower spectral resolution for 
clarity. 

spectropolarimeter on CFHT on August 1** and 
S'''* of 2006 and with HiVIS on November 7*'' 2006. 
We also have a back-to-back observation set of 
MWC361 and HD163296 with HiVIS on June 18- 
24th 2007 and with ESPaDOnS June 26th and 
27th. All data sets were taken in fair weather 
and reduced with the dedicated, automated Esprit 
data reduction package provided by CFHT (Do- 
nati et al. 1999). As a further check, we adapted 
the HiVIS reduction package to process the ES- 
PaDOnS data and found very similar results with 
both packages. 

The ESPaDOnS data shows that over a period 
of three nights in 2006 (Aug 01 to Aug 03) and 
again over 2 nights in 2007 (June 23rd and 24th), 
the spectropolarimetry from MWC361 is identical 
within the noise. The spectropolarimetry nearly 
a year apart matches very closely, though the in- 



tensity profiles change slightly. Our HiVIS obser- 
vations, while at lower spectral resolution, show a 
similar polarized spectrum. A side-by-side com- 
parison of the MWC361 data is shown in figure 
HOI Each data set from HiVIS has been rotated 
by an arbitrary angle to take out any telescope 
frame rotation and plane of polarization rotations. 
These by-eye rotations align the polarization spec- 
tra quite well, showing that HiVIS can reproduce 
the ESPaDOnS measurements quite well. 

The ESPaDOnS observations of HD 163296 
show a large amplitude polarization change, al- 
most 1% peak-to-peak, near the absorptive com- 
ponent of the emission line. Observations with 
HiVIS also show a similar effect at these same 
wavelengths, but the magnitude of the change in 
the HiVIS data is smaller. Figure |41] shows the 
side-by-side comparsions of HiVIS and ESPaDOnS 
data. The HiVIS data has also been de-rotated 
but the fits are not as good for many of the data 
sets. This target had a different and wider range 
of pointings and could have the polarization am- 
plitude modulated by the telescopes polarization 
properties. Still, one can clearly see the proper 
qu shape and even similar values in the HiVIS 
observations. 
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Fig. 42.— The spectropolarimerty of MWC361 
from CFHT on Aug 01 and 03 2006 plotted as 
degree-of-polarization and position angle. The 
continuum polarization has been shifted to match 
the Mottram et al. 2007 polarization of 0.82% at 
95.6° or q-u continuum of (-0.80%, -0.16%) before 
calculating these polarization spectra. 

We can also compare our observations to those 
of MWC361 taken with the Wilham Herschel Tele- 
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scope (WHT) using the ISIS spectropolarimeter as 
presented in the Vink et al. 2002, 2005b and Mot- 
tram et al. 2007 papers. The WHT is an alt-az 
telescope and ISIS is a long-slit spectropolarime- 
ter. The slit is aligned on the sky by physically 
rotating the entire spectrograph. The polarizing 
optics are almost identical to the HiVIS design 
- a Savart plate just after the slit creates the to 
orthogonal polarizations with a rotating quarter- 
and half-wave retarder mounted in front of the 
slit. A polarization observation is four sets of four- 
wave-plate rotations. The slit however has only 
the primary and secondary mirrors upstream, and 
those near-normal reflections do not induce sig- 
nificant instrumental polarization (Harries et al. 
1996). However, the waveplate does produce a 
very significant 0.2% amplitude ripple that must 
be subtracted from every high-resolution data set. 




Fig. 43.— The qu- loops from MWC361 from 
CFHT on Aug 01 and 03 2006 shown with the 
Mottram et al. continuum of -0.80%, -0.16%. 
Stokes q decreases broadly over the line and u pro- 
vides the small vertical deviations at the qu-loop's 
maximum amplitude. 

These WHT observations had a resolution of 
35km/s or R~8500, whereas the 1.5" slit on the 
AEOS spectrograph has a resolution of R^13000 
in the typical polarimetric mode. All data sets 
are then binned by flux to a constant polarimetric 
error. Vink et al. 2002, hereafter V02, present 



data from December 1999, Vink et al. 2005b, 
hereafter V05, presents data from December 2001 
and Mottram et al. (2007) present data from late 
September 2004. All three data sets show a simi- 
lar qu-loop and polarization spectrum, though the 
Vink et al. 2002 data has much higher noise. The 
continuum polarization was reported as 0.80% to 
0.82% at a PA of 94° to 96° in V02 and V05 re- 
spectively. We can take our CFHT and HiVIS 
data and scale them to these continuum values 
(q=-0.8% u=-0.15%) and calculate the degree of 
polarization and PA. We find a very similar po- 
larization spectrum, though with a significantly 
different PA spectrum, shown in figure [42l This 
changing PA is caused almost entirely by the u 
spectrum changing sharply in the center of the 
line while q is positive over the whole broad re- 
gion around the line. In the V02 and V05 mea- 
surements, the PA change is on the red side of 
the emission line, with a roughly 20° change. Our 
measurements now show a double-peaked change, 
still on the red side of the emission line with a 
similar wavelength range, but now with a two-fold 
change: a 15° increase in PA followed by a 15° 
decrease. The significance of this is best seen in 
the qu-loops, shown in figure 1151 Our plot shows 
the u changes as both positive and negative with 
extreme values near -0.4% and 0.1% whereas the 
V05 qu-loops show only a significant decrease in u 
with a barely noticeable increase. 

This illustrates the type of conclusions one can 
draw with HiVIS data. MWC361 has a high dec- 
lination, -|-68, and the range of pointings for all 
observations were quite small. The star is always 
low and in the north, from altitudes of 30-42 and 
azimuths of 340 to 20. HiVIS was quite successful 
at reproducing the spectropolarimetric measure- 
ments between ESPaDOnS and ISIS. HD163296 
however has a declination of -21 and spans a much 
greater range of pointings - altitudes from 30-48 
and azimuths of 132-226. The signature is still de- 
tected at the correct wavelengths, but at a reduced 
amplitude. 

Another comparison between the William- 
Wehlau spectropolarimter or WWS and HiVIS 
is possible (Eversberg et al. 1998). The WWS 
has a design very similar to ESPaDOnS, but us- 
ing two achromatic quarter-wave plates instead 
of two half-wave and one quarter-wave Fresnel 
rhombs. Eversberg et al. 1998 report commis- 
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Fig. 44. — Twenty-five q, u, and intensity spec- 
tra of 7 Cassiopeia taken on many nights from 
September 2006 to September 2007. The plane 
of polarization has been rotated for each spectra 
to maximize q in the core of the emission line to 
provide a relative reference frame. The average 
of all the polarization spectra is illustrated as the 
thick solid line showing a broad 0.2% polarization 
change across the line. 



sioning observations of 7 Cas on the 1.6m Obser- 
vatoire du mont Megantic (OMM) telescope that 
reproduce the 0.2% polarization change originally 
seen in a Pockels-cell polarimeter measurement 
with the University of Western Ontario 1.2m tele- 
scope (Poeckert & Marlborough 1977). We see a 
similar smooth and broad signature in 7 Cas as 
well, shown in figure [331 The observations have all 
been rotated by an arbitrary angle to maximize 
the -|-q spectrum in the core of the emission line. 
With the simple, broad form of the polarization 
signature, this has the effect of producing an effec- 
tive de-rotation to a common polarimetric frame. 
The resulting polarization spectra are averaged to 
produce our best estimate of the average polar- 
ization signature, shown as the thick line in figure 
l44l This average signature shows an amplitude of 
nearly 0.2% with a very smooth profile. 

In summary, HiVIS is capable of reproduc- 
ing spectropolarimetric measurements made with 



many other instruments. However, there are some 
complications that arise from the polarization 
properties of the coude path that cause some spec- 
tropolarimetric measurements to be significantly 
reduced in magnitude and/or rotated. 

6. Conclusions 

We have presented our dedicated IDL-based re- 
duction package and the polarization properties of 
the HiVIS spectropolarimeter. This package com- 
pares very favorably with other dedicated reduc- 
tion packages. The Esprit package on ESPaDOnS 
and our scripts give very similar results when ap- 
plied to ESPaDOnS data. We have found that 
the telescope induces 0.5% to 3.0% continuum po- 
larization in an unpolarized source with a strong 
dependence on secondary focus. The telescope 
can also strongly depolarize a source or rotate 
the plane of polarization as was shown by highly- 
polarized twilight observations. All of these are 
functions of altitude, azimuth, and wavelength. 
Optical ray trace models with Zeemax were pre- 
sented to illustrate the pointing dependence of the 
telescope polarization properties, but were ulti- 
mately shown to be quantitatively inaccurate due 
to the very common problem of oxide coatings. 

While these telescope effects complicate the 
analysis of data taken with this instrument, use- 
ful information about the amplitude and wave- 
length dependence of the polarization across a 
spectral line can still be obtained. The perfor- 
mance of the instrument was compared to other 
spectropolarimeters, particularly ESPaDOnS on 
CFHT, ISIS on WHT, and WWS on OMM. HiVIS 
found to reproduce the spectropolarimetric signa- 
tures, sometimes strikingly well given the possible 
complications. This instrument will be able to 
make sensitive spectropolarimetric observations, 
despite some complications in the interpretation 
of the polarization after so many oblique reflec- 
tions. 
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